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ABSTRACT 



An interactive computerized experimental procedure is 
presented for obtaining the elastic-plastic fracture 
toughness <Jj c ) for a material. The process employs a bend 
type, single specimen approach by using the measured 
unloading crack opening displacement (COD) compliance to 
obtain the crack length. Data are taken, manipulated and 
plotted during the test to provide the operator with the 
necessary information to control the testing machine. The 
program automatically plots the J versus crack extension 
curve and determines Jj c upon completion of the test. A 
second program is presented to assist the operator in 
fatigue precracking the specimen. A third program is used 
to recall, list and plot the data that were stored on 
diskette during the test. 
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INTRODUCTION 



A. GENERAL 

The toughness of a Material is a Measure of the energy 
per unit volume required to cause fracture. It can be 
quantified by simply taking the area under the true stress 
strain curve from a siMple tension test [Ref . 13 . The 
practical usefulness of this aeasureaent is liaited however, 
because it provides no inforaation on the material's ability 
to resist crack propagation. Since all real structures 
contain notches or cracks of some sort, by design ox- 
otherwise, which act as stress concentrators, another 
measurement of toughness was developed called fracture 
toughness. 

The fracture toughness <K C ) of a material is a Measure 
of the resistance of the Material to fracture in the 
presence of a defect. K is defined as the plane strain 
fracture toughness and is used to quantify a Material's 
resistance to crack propagation [Ref. 21 . The American 
Society for Testing and Materials (ASTM) developed 
Standard E399-83 [Ref. 3] which establishes precise 
details for accepted testing procedures. The standard 
also provides rigid criteria for the specimen size and 
geometry to ensure plane strain conditions exist at the 
tip of the crack. 
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The J-integral method, developed by Rice in 1968, 
provides a means o£ including plastic deformation at the 
crack tip [Ref. 43. J is a measure of the intensity of the 
entire elastic-plastic stress-strain field that surrounds 
the crack tip. Jic is defined as the value of J at the 
onset of cracking. J-integral testing has gained acceptance 
as a valuable means of evaluating a metallic material's 
resistance to crack initiation as demonstrated by 

establishment of ASTM Standard E813-81 CRef. 53. The 
measured value of Ji c has practical applications in 

measuring the effects of metallurgical variables, heat 
treatments, and weldments and is being used by researchers 
today to make such quantitative distinctions [Ref. 5-73. 

B. Jic VERSUS Kic TESTING 

The procedures established by ASTM Standard E399-83 for 
obtaining Kic impose minimum specimen size requirements 
that cause difficulty when testing materials that possess 
high toughness. 'The specimen thickness must be sufficient, 
to cause plane strain conditions to exist [Ref. 33. 
Hertzberg [Ref. 23 presents the empirical equation 
developed by Brown and Srawley [Ref. 83: 

t and a > 2.5 * <Ki c /Ys3 

where 

t * specimen thickness 



a 3 crack length 



Kic 3 plane strain fracture toughneaa 
Ya 3 yield etrength 

A aimilar analyaia for J-integral teat apecimena reveala 
that : 



where 



b or B > 25 » (Ji c /Fa> 



b 3 initial uncracked ligament 

B 3 specimen thickness 

Jlc 3 anticipated value if Ji c 

recheck when teat ia completed 
Fa 3 flow atreaa [ ( Ya+Uta) /23 



Hertzberg takea thia reault one atep further by ahowing that 
Jlc equala <Kic>^/E (where E equala the material's elaatic 
modulua) for Jxc« By dividing the two equationa one aeea 
that: 

(Kxc specimen size) > 20 * (Jlc specimen size) 
thus, because the J-integral allowa plaatic deformation, one 
of its moat significant advantages ia that a satisfactory 
specimen can be obtained at one- twentieth (1/20) of the 
thickness of a Kic specimen. Thia not only meana smaller 
specimens may be used but, alao, smaller testing machines 
are required and therefore coats are reduced. CRef. 23 

Although Jic and Kic provide a different measure of 
toughness, Jic can be used to estimate Kx c . ASTM Standard 
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E813-S1 CRef. 53 provides United guidelines on the 
relationship between the two. 

C. PURPOSE AND SCOPE OF RESEARCH 

The purpose of this thesis is to develop the procedures 
required to obtain Ji c fron J-integral test results using 
equipment currently available within the Materials Science 
group of the Mechanical Engineering Department at the Naval 
Postgraduate School (NPS) . A list of the major equipment 
follows: 

a. Series 810 MTS machine 

b. HP-9826 computer with internal disk drive 
and 320K random access memory 

c. HP-3497A data acquisition system 

d. HP-3437A voltmeter 

e. HP-2671G thermal printer 

f. HP-7225B graphics plotter 

Figure 1 shows the system layout. > 

The teat method developed is to comply with ASTM 
Standard E813-81 CRef. S3, be easy for others to use, and 
require a minimum of operator interaction. It is 
specifically designed for use with three point bend 
specimens using a single vice multiple specimen approach. 
The method developed and presented here covers the entire 
test, including specimen preparation, precracking, testing. 



and validation of results 



II. MATURE OF THE PROBLEM 



A. J-INTEGRAL 

The J- integral is a mathematical expression uaed to 
characterize the near tip atreaa-atrain environment o£ 
cracked elaatic-plaatic bodies [Ref. 9] . It ia a path 
independent line integral that surrounds the crack tip for 
non-linear plastic metallic materials CRef .53 . The equation 
delveloped by Rice is CRef. 4] : 



W * loading work per unit volume, or for elastic 
bodies, strain energy 

(~* 9 the path of the integral, which encloses the 
crack tip 

ds 9 increment of the contour path 
T 9 outward traction vector on da 
u 9 displacement vector at da 
x,y,z 9 rectangular coordinates 

T‘(£u/Ax)ds 9 the rate of work input from the stress 



Figure 2 provides a pictorial representation of this 
equation. 

The mathematical expresaion developed by Rice CRef. 41 
has been manipulated and reduced to simpler forma for 




where 



field into the area enclosed by f* . 
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specific specimen geometries. Rice, Paris, and Merkel [Ref. 
10] reduced the J-integral expression for a three point bend 
specimen to: 

2 * A 

J 

b * B 



where 

J = resistance to crack initiation 
A = area under load versus load line 
displacement curve 
b s uncracked ligament 
B * specimen thickness 

This thesis is not intended to redevelop the J-integral but 
rather apply it, therefore only the result is presented 
here . 



More recent work by Ernst, Paris, and Landes CRef. 11] 
corrects the calculation of J due to crack growth by means 
of the following equation: 

J(i+1> * < Ji* <n/b) i*Ai , i+i ) * < 1- < Y/b) * (ai+i-ai ) ) 

where 

J<i+1> = current value of J 
Jj. = last J computed 

n * 1 for three point bend type specimens 
Ai,j. + i s area enclosed by the actual test 
record and lines of constant 
displacement dj. and dj. + i 
Y s 1 for three point bend type specimens 
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b 3 current uncracked ligament 
W = specimen width 

= last crack length 
ai-t-1 3 current crack length 

Correcting for crack growth during the teat is necessary 
when using a single specimen approach as used by this 
method . 

The blunting line defined as: 

J » 2 » Fa * delta a 



where 



Fa s flow stress 

delta_ a 3 change in crack length 
approximates apparent crack advance due to crack tip 

blunting in the absence of slow-stable crack tearing. It is 
plotted on the J versus crack extension plot. A linear 
regression fit through the applicable data points is then 
plotted. The intersection of the two lines is by definition 
Jjc (see Figure 3) . CRef . 51 

B. PREVIOUS WORK 

The early experimental methods used to obtain Jjc 
required the preparation of several specimens. Each 

specimen received a different load. The specimens were then 
heat tinted (to identify crack length) and broken apart. 
Each specimen provided one data point on the J versus crack 
extension curve, commonly referred to as the J-R curve CRef. 
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9 and 11) . Checks were then performed on the validity of 
the data. A minimum of four valid data points were required 
to obtain a linear least squares fit and hence Ji c [Ref. 
5,10 and 12] . 

Researchers over the last five years have developed a 
much improved test method whereby a single specimen can be 
used. Today's computerized data acquisition equipment and 
data manipulation are ideally suited for the task. By 
relating the compliance (how much a specimen bends under a 
given load) to the crack length, all the data required to 
obtain the J-R curve can be collected by using a single 
specimen [Ref . 11-16] . Clarke [Ref. 16] states that prior to 
the use of computers the time required for data reduction to 
obtain the J-R curve was four to five hours. Using 
computerized, single specimen, unloading compliance methods 
it is now possible to to obtain the J-R curve in a matter of 
twenty to thirty minutes after the test starts. The obvious 
savings in specimen preparation and data reduction have made 
the J-integral a valuable means of analysis where no 
practical means previously existed. 

C. TEST CRITERIA 

The test method developed in this thesis is specifically 
designed for use on a three point bend specimen. It 
utilizes the single specimen approach by relating the crack 
opening displacement (COD) compliance to the crack length. 
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The following equations are currently used by DTNSRDC. For 
bend specimens with a span to width (S/W) ratio of 4.0, Dr. 
J. Joyce CRef. 17] of the United States Naval Academy 
determined: 

a/W * 0.998265 - 3.81662 * U - 1.80596 * U2 ♦ 32.31038 » U3 
- 44.15665 • U4 - 52.67876 » US 

where 

a * crack length 
W * specimen width 



U = 



e 



Be 



W 



compl 



^ 1/2 



where 

B e 3 effective thickness = Bmax - (Bmax - Bmin)2/Bmax 
E * material's elastic modulus 
S 3 span between supports 
W 3 specimen width 

compl 3 compliance 3 1/slope of the load versus COD line 
for the specific unloading. 

The reason for using the single specimen approach has 
already been presented. The reasons a bend type specimen 
and COD type unloading compliance are used are: 

a. Three point bend specimen geometry is better suited 
than compact tension specimens for testing weldments. 

b. Three point bend specimen geometry is also easier to 
machine than compact-tension specimens. 
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c. COO voltage readings are less susceptible to noise 
than linear variable differential transformer (LVDT) 
voltages used to measure load line displacement . 

d. The COO compliance equation affords 0.1H fitting 
accuracy and allows comparison of results with 
DTNSRDC test data CRef. 173. 
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III. TEST METHOD 



A. OVERVIEW 

A J-integral teat can be divided into three distinct 
steps, 1) pre-test, 2) test, and 3) post-test. Each o£ 
these are described in greater detail in this chapter and 
presented as a step by step procedure in Appendix A. 

After the specimen is machined it must be precracked 
using cyclic fatigue loading to minimize residual stresses 
and create a sharp crack tip CRef . 53 . A computer program 
was developed to assist the operator in accomplishing this. 
The program computes the loading parameters for precracking 
in accordance with ASTM Standards E399-83 and ASTM E813-81 
CRef. 3 and 53 . The operator uses this information to run 
the MTS machine. 

During the J-integral test a second computer program 
collects and records the load, COD, and load point 
displacement (deflection) . These data are used by the 
program to determine J and the amount of crack extension for 
each operator controlled unloading cycle. An unload cycle 
consists of decreasing the load by 10* of the expected 
maximum load and reloading to the original load point. The 
program provides the operator with a real-time graphical 
representation of the any two parameters (load versus COD or 
deflection) during the test. It will also plot the J-R 
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curve and compute Ji c upon completion of taking all the teat 
data. Data are stored on diskette to serve as a record and 
allow future recall listing or plotting. 

The post-test process includes heat tinting the specimen 
to aid in visual verification of crack extension upon 
completion of the test as required by ASTM Standard E813-81 
CRef . 53 . 

B. SPECIMEN PREPARATION 

Figure 4 shows specifically the dimensions used in 
preparing eight specimens from two different one inch thick 
welded plates of HY-80. One plate was preheated while the 
other was not. It was necessary to reduce the width to 
approximately 0.95 inches to remove the distortion (bow) 
created by the weld. Four specimens were prepared from 
each plate, the machined notch was placed in the weld metal 
for two and in the base metal for the other two. The clip 
gage (used to measure COD) is attached to the top of the 
specimen by means of two cut razor blades. The blades are 
held down by small straps which allow adjustment of the 
initial gap. The initial gap is set to 0.475 inches, the 
calibrated zero setting for the clip gage, by means of a 
machined gsge tool that rests in the specimen's machined 
notch. This method is also used by DTNSRDC. 

Two additional holes are drilled and tapped for 
retaining a small bracket on the face of the specimen. The 



bracket is used to provide a sensing point for a linear 
variable differential transformer (LVDT) . The LVDT signal 
is used to calculate the area under the load versus 
displacement curve. At present the external LVDT is unable 
to control the MTS machine so the current method uses the 
MTS' internal LVDT to measure deflection via actuator 
movement . 

Clarke [Ref. 16] states that side grooves improve the 
accuracy of calculating the crack length when using 
unloading compliance procedures. The side grooves reduce 
plastic deformation on the surface and aid in propagating a 
straighter crack front. DTNSRDC has found 20X side 
grooving produces the best results, hence the same was used 
for this work. 

After the specimen is prepared a test record sheet (see 
Appendix B) is completed. Notch and side groove depths can 
be measured using a dial indicator with a needle tip. The 
record provides the operator with ready access to the 
information necessary to run the precracking and test 
programs . 

C. PRECRACKING 

Before the J-integral test can be performed the 
specimen must be precracked to establish a sharp crack tip. 
ASTM Standards ES13-81 and E399-83 establish parameters that 
must be met during precracking to prevent over-stressing the 
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specimen or propagating the crack through the entire width. 
Precracking is accomplished with the aid of a program called 
"PRECRACK" which will compute loading requirements based 
upon these paramenters. 

The program, written in HP-Basic, affords flexibilty by 
allowing the use of different bend specimen sizes and 
materials. It is easy to use due to its interactive design 
and easy to modify should such a need arise. A listing of 
"PRECRACK" is provided in Appendix D. 

The operator is requested by the program to enter 
specific specimen data. The program will compute the initial 
maximum and minimum loads and change in the associated 
stress intensity factor (delta K) . If this delta K exceeds 
40 ksiftinl/^ then new maximum and minimum loads are 
computed. The operator is responsible for ensuring the MTS 
machine is operating within the prescribed limits. The 
operator can obtain a hard copy of the entered and computed 
data for record purposes (see Appendix G) . 

The MTS machine produces cyclic loading by means of an 
inverted (compressive load) haver sine wave. 

The operator must periodically interrupt the cyclic 
loading and measure the current crack length by manually 
loading and unloading the specimen (determine the specimen's 
compliance) . The program computes the new crack length and 
new loading limits. During the last 50X of precracking: 

Pmax < . 4 Pl 
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where 



Pmax 3 maximum load 

PL ■ t <4/3> * <B * b2 * Fa / <2 » Hs))3 

where 

B 3 net thickneaa (Bmin if specimen has side grooves) 
b 3 current uncracked ligament CW-a) 

Fa 3 flow stress ( (Yakuts) /2J 
Hs = half span 

and Delta_K/E < 0.005 mml^2 <0.001 inl^2) whichever, 

produces the least load CRef . 5] . The program will 

automatically begin using these parameters when the most 
recently calculated crack length exceeds this 50* of the 
required extension to reach a/W 3 0.6. This value is used 
because it is the optimum crack length with which to begin 
the J-integral test [Ref. 51 . 

D. J-INTEGRAL TEST 

The J- integral test is performed by slowly applying a 
load to the specimen and periodically unloading it by 10* of 
the expected maximum load. The MTS machine is set up in 
atroke control with a very amall inverted ramp function 
resulting in a loading rate of 0.005 inch/minute. The 
operator must again interrupt the MTS program and manually 
unload the specimen. Each unloading cycle (unload and 
reload) provides the necessary information to compute the 
current crack length. The crack length is then used to 
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determine the value of J up to that point. The computer 
program "J_INTGL M provides the operator with pertinent 
information during the teat. Once the test is completed, 
the operator can obtain the J-R curve and value of Jj c (see 
Figure 4 and Appendix H> . The program will use only data 
points between the offsets, (blunting line ♦ 0.15mm) to 
(blunting line ♦ 1.5mm) as required ASTM Standard E813-81 
when computing the linear regression line. A minimum of 
four data points must be within the offsets and one must lie 
near (within 1/3 the overall range of delta_a> the blunting 
line. Specific acceptability criteria for data is provided 
in ASTM Standard E813-81 CRef. 51. 

The program M J_INTGL" is the foundation of the test 
method, the single specimen approach would not be possible 
without it. The program collects, manipulates, displays and 
records all pertinent data during the test. A flow chart of 
"J_INTGL‘* is provided in appendix C. A listing of M J_INTGL” 
is provided in Appendix E. 

Deflection measurements are presently taken using the 
LVDT installed in the MTS machine. A correction is 
therefore required to account for the system's compliance 
while loading a specimen. The correction was obtained by 
loading the system with a two inch thick piece of HY130 and 
recording the stroke reading for various loads. This 
process was checked three times with consistent (within 5X) 
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results. A plot of the dats and the 4th order ploynomial 
curve fit is included as Figure 5. 

E. POST TEST VERIFICATION 

The results must be verified upon completion of the 
test. ASTM Standard E813-S1 presents the following criteria 
that must be met before the test can be considered valid: 

1. Total crack extension (visually measured) must be 
within 15x of that recorded during the teat. 

2. bi and Bain > 25 • Jj c / Fs (all terms are as 
defined earlier except bi 3 post-fatigue uncracked 
ligament) . 

3. The slope of the linear regression fit to the data 
must be less than flow stress (dJ/da < Fs) . 
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IV 



RESULTS 



A. TEST RESULTS 

Specimen HY80-5B, prepared from the base metal of the 
preheated welded plate, was precracked, and tested resulting 
in a Jlc value of 142 kJ/m^. Figure 8 is a tracing of a 
photograph of the final fracture surface after heat tinting. 
All post teat checks were completed and yielded satisfactory 
results. The Ji c value is in good agreement with previous 
tests of HY80 plate as conducted by DTNSRDC resulting in an 
average Jic of 144 kJ/m^ [Ref. 6] . 

Plots of J-R curve, load versus COD, load versus 
deflection and the are included as Figures 5, S and 7, 
respectively. The load versus deflection curve was obtained 
using the program * , J_DATA", thus proving that the data is 
stored on disk, and is able to be recalled. 

Data was taken at a nearly constant rate of 2.5 data 
points per second (150 per minute) while the program was 
executing . 

Test reults were verified following the guidelines 
provided by ASTM Standard E813-81 . The fracture surface was 
enlarged by taking a photograph, a tracing of which is 
included as Figure 8. The observed crack length was 14. 4* 
longer than the last recorded crack length. This is within 
the acceptable limit (recall ASTM requires < 1550 but 
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because the clip gage reached its limit detector and caused 
the test MTS to stop between unloadings, it is possible some 
crack extension occurred after the last unloading cycle. 
This would help to explain the difference. 

Additional information on this specific test is provided 
by being used in the example calculations shown in Appendix 
A. The exemplary computer output for each of the programs 
is also from specimen HY80-5B. (See Appendixes G, H and I). 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The method developed and presented by this work produces 
accurate and reliable results that meet the specifications 
established by ASTM Standard E813-81. Precracking specimen 
HY80-5B took approximately 2 hours and the J-test required 
approximately 2.5 hours including several pauses to manually 
verify several calculations. The time required for the 
processes will most likely decrease as operator experience 
and confidence increases. 

All three computer programs have run successfully and 
perform as intended. 

B. RECOMMENDATIONS 

As with all first versions of computer programs as 
extensive as "J_INTGL", changes and improvements are to be 
expected. The following programming changes may be of value 
to future operators: 

1. Display pertinent information (i.e. current crack 
length in both English and metric units (currently 
only metric is provided) . 

2. Provide a means of storing the data arrays during the 
test without being required to rerun the entire 
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program (re-entering specimen data 



conversion data 



etc . ) 

3. Provide a means of eliminating a specific J-R data 
point from being used in obtaining the linear 
regression fit. This may be programmed to be 
accomplished automatically by applying Chauvenet's 
Criterion to the data. 

Improvements in the current means of conducting the test 
can also be made, such as making the MTS control both 
loading and unloading using the same rate, requiring the 
operator to do no more than press a button or throw a 
switch. Lastly, pursue the use of the external LVDT, this 

will eliminate the need for correcting the deflection 
reading and speed up the rate of data acquisition. 
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MTS 810 




Figure 1. System Layout 
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J-Integral 




Figure 2. Pictorial Representation of J-Integral 
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PLOT of J vs CRACK EXTENSION for HY80-5B 
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Figure 3. J Versus Crack Extension for HY80-5B 
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Figure 4. Specimen Drawing 
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SYSTEM COMPLIANCE 
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Figure 5. Plot of System's Compliance 
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Figure 6. Plot of Load Versus 1 COD for HY80-5B 
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Figure 7. Plot of Load Versus Deflection for HY80-5B 
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Figure 8. Tracing of HY80-5B Fracture Surface 
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APPENDIX A 



TESTING METHODS 



A. INTRODUCTION 

The steps required to perform a J-integral test using a 
single specimen approach can be divided into three major 
categories, pre-test, test, and post test. Each of these 
categories may be further divided as shown below. 

1. Pre-test 

a. Machine specimen 

b. Precrack specimen by fatigue loading 

2. Test 

a. Run J-integral test 

3. Post- test 

a. Heat tint specimen and break apart 

b. Measure crack extension and verify results 



Each process 


is described in 


greater 


detail in 


thia 


appendix. 


A 


working knowledge 


of the 


equipment 


is 


presupposed 


and 


awareness of all 


safety 


precautions 


ia 



required before operating the MTS. 

B. PRE-TEST 
a. Machine Specimen 

ASTM Standard E813-81 provides very specific guidance on 
specimen dimensions. All dimensions are related to the 
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width (W> . Fatigue precracking ia enhanced if the radiua at 
the tip of the notch ia kept to approximately 0.003 inchea. 
Thia waa accomplished by machining with a apecial cutter 
which waa uaed to make only the final paaaea in the notch. 
Figure 2 can be conaulted for the actual dimenaiona uaed to 
machine the apecimena uaed in thia work. 

b. Precracking 

The following atep by atep guide ia provided to aaaist 
operators . 

1. Before turning on the MTS machine ensure the proper 
range card ia inatalled in the atrain aignal conditioner for 
the clip gage to be used. 

2. Turn on power to all equipment with the computer being 
switched on laat. 

3. Load and run "PRECRACK" program on the HP-9826 computer. 

4. Select range acalea on NTS; 

example: 10 v * 4000 lba 20X of full load (range #3) 

10 v = 0.02 in 10X of full clip gage (range #4) 

10 v * 2.5 in SOX of full stroke (range #2) 

5. Preaa kO-key on computer and enter specimen data as 
requested . 

6. Preaa k2-key and enter conversion data aa requested . 

7. Preaa k5-key to obtain a table of COD readings for 
various loada and crack lengths. 

8. Set READOUT knoba on the MTS Controller to be: 
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top = load 
Biddle = atrain 
bottom * stroke 

Set all to <♦) polarity. See Figure Al. Thia enaurea 

channela of the MTS machine are properly aligned with the 
data acquiaition ayatem. 

9. Set MTS limit detector awitchea baaed on the information 
juat generated by the program. Limit awitch aettinga are 
determined by the aimple formula: 

Limit awitch * L • 10 / FSR 



where 

L = deaired limit in Engliah unita 
FSR = full acale range of current operating 
being uaed in Engliah unita 
example "PRECRACK** determined: 

Pmax = 1668 lba aet load limit to 1600 lba 
for apecimen HY80-5B 

load limita: 

upper limit amall poaitive load for example 40 lba 

<ahould never be poaitive) 

♦40 lba * 10 v / 4000 lba = ♦O . 1 v 
therefore aet upper limit knob to +0.1 v 

lower limit 1600 lba maximum load 

-1600 lba » 10 v / 4000 lba = -4.0 v 
therefore aet lower limit knob to -4.0 v 
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strain limits: 



upper limit 0.01 in beyond maximum expected during 

precracking 

•♦■0.01 in * 10 v / 0.02 in = ♦5.0 v 

lower limit some small value for example -0.1 v 

(should never be negative) 

stroke limits: 

upper limit 1.5 in to allow clearance for installing 

specimen and yet stop the MTS if this 
limit is reached 

♦1.5 in * 10 v / 2.5 in = ♦S.O v 

lower limit 0.25 in of compressive actuator movement 

is more than sufficient for precracking 
-0.25 in * 10 v / 2.5 in = -1.0 v 

10. Adjust the upper roller supports as required to obtain a 
span-to-width ratio <S/W) of 4.0. 

11. Install specimen (with the clip gage already installed) 
using alignment collars to properly center. 

12. Zero the load and strain voltages. 

13. Reset and clear all limit switches (all green lights on 
controller panel should be showing) . 

14. Using load control energize hydraulics and apply a low 
(<10O lbs) compressive load to the specimen by SLOWLY 
turning the SETPOINT knob counter-clockwise. 
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15. Slowly decraae the coRpreaaive load on the specimen to 
between 5 and 10 lbs by turning the SETPOINT knob clockwise 
and then zero the stroke voltage reading. 

YOU ARE NOW READY TO OBTAIN AN UNLOADING COD COMPLIANCE 
MEASUREMENT OF THE CRACK LENGTH 

16. Press k4-key on the computer and proceed with taking 
data . 

17. Slowly, and as smoothly as possible increase the load to 
approximately 500 lbs and return to approximately 20 lbs 
ending the loading cycle with a alight (10 lbs) increase to 
trip the computer out of the data-taking loop. 

The computer will display all pertinent crack length 
information and ask if you want another measurement. You 
may repeat this process as often as you like, just be sure 
the maximun load is lower than that required to extend the 
crack . 

18. Make necessary adjustments to the MTS such that the 
strain limit switch will trip the MTS when the crack extends 
to SON of the desired extension. This is determined by 
interpolating results from the table obtained earlier. (See 
Step 7) . 

NOTE: As the crack propagates the acceptable maximum 
load decreases. For the specimen, HY80-5B, used to test 
this program a Pmax = 1000 lbs was suitable to initiate the 
crack and propagate it to SON of the desired extension 
(a = 0.506 in). 
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Example : 

SOX crack extension * crack length (a) of 0.506 in 
from table: 

COD * 1.7 volts at 1000 lbs and a 3 0.506 in 
Therefore set UPPER STRAIN LIMIT SWITCH to +1.7 v. 

YOU ARE NOW READY TO BEGIN CYCLIC LOADING 

19. Setup MTS: 

\ 

DISPLAY: 

a. Setup to read maximum and minimum load in lbs. 
FUNCTION GENERATOR: 

a. CONTROL MODE: Remote 

b. RATE l: 10 Hz 

c. RATE 2: N/A 

d. HOLD BUTTON: In 

e. BREAKPOINT: N/A 

f. HAVER SINE: In 

g. INVERT: In 

CONTROLLER : 

a. SPAN l: 0.0 

b. SPAN 2: 0.0 

COUNTER: 

a. COUNTER INPUT: Oscillator 

b. COUNTER MULTIPLIER: X100 

c. Press RESET 

20. Slowly increase the SETPOINT to Pmin. 

21 . Press the RUN button on the CONTROLLER and then release 
the HOLD button of the FUNCTION GENERATOR. 

22. Slowly increase the span of the cyclic loading using 
SPAN 1 knob to the desired Pmax. 
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The MTS will continue cyclic loading until it ia 
manually stopped or a limit detector ia tripped. Most 

likely the COD limit detector will trip at which time the 
machine will will stop and unload the specimen. 

23. A new crack length measurement can now be obtained 
using the same manual process as described in Step 16 and 
Step 17. 

24. Readjust limit switches and resume fatigue loading as 
done in Steps 18 through 22. 

25. Continue this process until the crack reaches the 
desired length. 

C. J- INTEGRAL TEST 

The following steps are required to run a J-integral 
test. They are very similar to those used in precracking. 
The major change ia that a ramp function under stroke 
control ia now used to control the MTS machine. 

1. Before turning on the MTS machine ensure the proper 
range card is installed in the strain signal conditioner for 
the clip gage to be used. 

2. Turn on power to all equipment with the computer being 
switched on last. 

3. Load and run "J_INTGL” program on the HP-9826 computer. 

4. Select range scales on MTS; 

example: 10 v = 4000 lba 20X of full load (range #3) 
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10 v - 0.02 in 10X of full clipgage (range #4) 

10 v * 0.5 in 10X of full stroke (range # 2) 

5. Press kO-key on computer and enter specimen data as 
requested . 

6. Press kl-key and enter conversion data as requested. 

7. Press k3-key and follow instructions for creating a data 
file on an initialized data diskette. See HP-Mannual on how 
to initialize a diskette if necessary. 

S. Press k3-key to setup the plotter. 

9. Press k5-key and obtain the maximum load you may apply 
while measuring aO. 

YOU ARE NOW READY TO SETUP THE MTS MACHINE 

10. Set READOUT knobs on the MTS Controller to be: 

top = load 
middle = strain 
bottom = stroke 

Set all to (♦) polarity. See Figure A1 . 

11. Set MTS limit detectors switches, 
example: 

Pmax = 895 lbs for HY80-5B (see Appendix H> 
load limits: 

upper limit small positive load (40 lbs) 

(should never be positive) 

•►40 lbs * 10 v / 4000 lbs = ►0.1 v 

lower limit 895 lbs compression maximum 

-895 lbs * 10 v / 4000 lbs = -2.125 v 
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strain limits: 



upper limit want to be able to read full scale 

on range #2 <0.1 in) 

+0.1in * 10 v / 0.1 in * +10 v 

lower limit a very small value set to -0.1 v 

C should never be negative) 

stroke limits: 

upper limit --- full scale <+10.0 v) on this most 

sensitive scale 

lower limit a maximum compression of 0.25 in is 

sufficient for the test and will stop 
the actuator before it can cause any 
serious damage to the specimen holder 
-0.25 in * 10 v / 0.5 in * -5.0 v 

12. Adjust the upper roller supports as required to obtain a 
span-to-width ratio of 4.0. 

13. Install the specimen <with the clip gage already 
installed) using the alignment collars to properly center. 

14. Zero the load and strain voltages. 

15. Reset and clear all limit switched <all green lights on 
controller panel should be showing) . 

16. Using load control energize hydraulics and apply a low 
<<100 lbs) compressive load to the specimen by SLOWLY 
turning the SETPOINT knob counter-clockwise. 

17. Decrease the load to between 5 and 10 lbs and zero the 
stroke voltage. 
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18. Adjust the MTS to the following settings: 

DISPLAY: 

a. Setup to read load In lbs and stroke in inches. 
FUNCTION GENERATOR: 

a. CONTROL MODE: Remote 

b. RATE l: 6000 <sec> 

as determined by: desire 0.005 in/min loading rate 

therefore ; 

0.5 in (full scale) / 0.005 in/mim • 60 aec/nin = 6000 sec 

c. RATE 2: N/A 

d. HOLD BUTTON: In 

e. BREAKPOINT: N/A 

f. RAMP BUTTON: In 

g. INVERT: IN 

CONTROLLER: 

a. SPAN l: 10.0 (full scale) 

b. SPAN 2: 0.0 

COUNTER: N/A 

YOU ARE NOW READY TO MEASURE THE INITIAL CRACK LENGTH 

19. The computer will begin taking data after the 'CONTINUE' 
key is pressed. The MTS will begin loading the specimen 
after the MTS' RUN switch is pressed and the HOLD button is 
released . 

20. When the desired maximum load (approximately 500 lbs) is 
reached; depress the HOLD button, switch to load control, 
and manually decrease the load (to approximately 20 lbs) by 

turning the SETPOINT knob clockwise. As with 
precracking you must end the unloading with a small increase 
in load (10 lbs from where you stop) to signal the computer 
you have completed the load/unload cycle. The computer will 
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plot, the date as you go, ensure you remain on the linear 
portion of the loading curve. 

21. The computer- will display all the pertinent crack length 
information. You may repeat this operation by switching 
back to stroke control and repeating Steps 19 and 20 as 
often as you like as long as you remain on the linear 
portion of the loading curve. 

22. After the computer has stored the desired initial crack 
length (a0> measurement on data diskette, you are ready to 
begin propagating the crack. 

23. Ensure you are in stroke control with near zero (<10 
lbs) load and the RUN light is on with the HOLD button 
depressed. 

24. Press the k7-key on the computer follow instructions and 
proceed to collecting data <data acquisition system will be 
continuosly switching channels and the plotter will be 
plotting), now release the HOLD button on the MTS. 

25. When you reach a point where you wish to unload the 
specimen; press the HOLD button, switch to load control, and 
slowly decrease the load. The computer will beep when you 
have unloaded the specimen by 10* of the expected maximum 
load. Switch back to stroke control and release the HOLD 
button to reload the specimen. 

26. When the load returns to the point from which the 
unloading began the computer will display and printout all 
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the information it haa determined with regarda to the 
current unloading cycle. 

27. Thia proceaa of awitching control modea (Stepa 25 and 
Step 26) are repeated for every unloading cycle until the 
you decide to atop by preaalng the k9-key on the computer or 
one of the limit detectors ahuta off the MTS machine. 
Remember, you muat have atleaat four data points between the 
offsets one of which lies near the blunting line on the J-R 
curve . 

28. The k9-key interrupts the computer program and tells it 
to stop taking data and write the data it haa collected on 
the data diskette. 

29. The computer will then display a table of the collected 
J versus delta_a and ask if you desire a plot. Simply follow 
the directions from the computer and you will obtain the J-R 
curve and the value of Jj c . 

D. POST-TEST 

The poat-teat includes verifying the final crack length, 
and validating the other parameters aa discusaed in ASTM 
Standard E813-81 auch aa: 

bi and Bmin > 25 * Jj c / Fa 

where 

bi * uncracked ligament after precracking 
Bmin » minimum specimen thickness 
Fa * flow stress 
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and that the slope of the linear regression fit to the data 
is leas than the flow stress (dJ/da < Fs) . This data is 
obtained directly from the "J_INTGL*‘ program output. 

Heat tinting the specimen, 10 minutes at 300C for most 
steels, aids in the visual identification of the amount of 
crack extension. After the specimen is removed from the 
oven and air cooled, place it in liquid nitrogen to assist, 
in obtaining a brittle fracture. After breaking the 
specimen apart, place the pieces in alcohol while returning 
to room temperature. This will prevent water vapor from 
condensing and freezing on the surface and minimize 
oxidation . 

A photographic enlargement of the surface can now be 
obtained and the crack length can easily be verified. The 
final crack extension must be within 15z of that recorded 
during the test . 
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Figure A1 . MTS Control Panel 
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APPENDIX B 



TEST RECORD SHEET 

THREE POINT BEND SPECIMEN 
TEST REPORT 

SPECIMEN I.D. MATERIAL ORIENTATION 

D«t» Operator T««t Taap 



T 

But ax 

i 



— T" 
Enin 



l 




W • A * 




b • 


L - 


Ha * 


Eaax * 




B*in • i 


Ya - 


UtA - 




1 

ft 

■ 


PolAAon' a ratio * 




Elaatlc 


RoduiuA • 



CoB««nU 



INTGL" FLOW CHART 



APPENDIX C 



J INTGL FLOW CHART 
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APPENDIX E 



PRECRACK LISTING 



in 
20 
30 
*0 
50 
60 
70 
80 
90 
100 
1 10 
120 
130 
1 40 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
250 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 



PRECRACKING PROGRAM FOR THREE °9INT 
BEND SPECIMEN. THIS PROGRAM INCORPORATES 
THE USE OF THE FOLLOWING EQUIPMENT: 

1. HP-9826 COMPUTER 

2. HP-3497A DATA ACQUISITION SYSTEM 

3. HP-2671G THERMAL PRINTER 

4. MTS MACHINE 

THE PROGRAM USES BOTH ENGLISH ANO METRIC 
UNITS. ALL METRIC VARIABLES ARE DENOTED 
WITH AN n» IN THEIR NAME le Wm-WIDTH IN 
METRIC UNITS <mm IN THIS CASE). 

THIS PROGRAM WAS WRITTEN BY: 

LCDR W. K. TRITCHLER. USN 
VERSION 1.0 DATE 1 1/33 



COM /J_intgl/ FI 

I 

I DIM Lodm<2000) ,Def lra<2000) «Codm<2000) these arrays will 
! required in the New a routine. 

! Lodm< I) - LOAD IN NEWTONS 

1 Deflm(I) - DEFLECTION IN nn 

! Codra< I) - CRACK OPENING DISPLACEMENT IN mm 



DIM Coef <6> .G<6) ,A0< 100) 
DUMP DEVICE IS 706 
Dvm-709 
CLEAR Dvm 



PRINT CHRS< 12) 
PRINT -• 

PRINT " 

PRINT " 

PRINT 
PRINT " 

PRINT 
PRINT " 



! CLEARS CRT 

PRECRACKING PROGRAM- 
VERSION 1.0" 

WRITTEN BY LCDR W.K. TRITCHLER” 
PLEASE SELECT A It -KEY*' 



Begin: ! 

ON KEY 0 LABEL "SPEC_DATA" GOSUB Spec_data 
ON KEY 2 LABEL ”CONV_DATA" GOSUB Conv_data 
ON KEY 4 LABEL “NEW_A" GOSUB New_a 
ON KEY 5 LABEL "a_TABLE” GOSUB A_table 



; 



» 



500 Spin: * WAITS FOR KEY INPUT 

510 GOTO Spin 
520 1 

530 ! * ** ***•••«*«««**»*-»*'•'•«'•*««***'•#***#**'•* * 

540 ! 



550 Spec_data: ! INPUTS REQUIRED SPECIMEN DATA 

560 PRINT CHR$< 12) ! CLEARS CRT 

570 PRINT "INPUT SPECIMEN IDENTIFICATION <ID>:" 

580 INPUT DS 

590 PRINT — 



be dimensioned 
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600 PRINT "J- INTEGRAL TEST Tor ,, :DS 
610 PRINT 

620 PRINT "INPUT SPECIMEN DATA USING ENGLISH UNITS (LBS. INCHES, etc.)" 

630 PRINT — 

640 PRINT "WHAT IS THE SPECIMEN'S WIDTH. NOTCH LENGTH, Bmax . Bmin (THICKNESS) 
?•• 

650 INPUT W.A.Bmax.Bnun 

660 PRINT — 

670 PRINT "WHAT IS THE SPECIMEN'S HALF SPAN AND POISSON'S RATIO?" 

680 INPUT Hs.Nu 

630 PRINT — 

700 PRINT "WHAT IS THE SPECIMEN'S ELASTIC MODULUS. HELD STRENGTH AND ULT 

IMATE STRENGTH?" 

710 INPUT E.Ys.Uts 

720 PRINT CHRS( 12) ! CLEARS CRT 

730 Fs” ( Ys+Uts) /2 

740 PRINT -USING ENGLISH UNITS ( in. ibs.psi)" 

750 PRINT "THE FOLLOWING HAS BEEN ENTERED for: ";DS 
760 PRINT "WIDTH NOTCH LENGTH -";A 

770 PRINT "Bmax *":Bmax:” Bmm »";Bmin 

780 PRINT "HALF SPAN POISSON'S RATIO *";Nu 

790 PRINT "YIELD STRENGTH »";Ys:” ULTIMATE TS -”;Uts 
800 PRINT "ELASTIC MODULUS »":E:" FLOW STRESS -";Fs 
810 PRINT — 

820 PRINT "THIS PROGRAM ALLOWS NO EASY RECOVERY, IF AMY DATA IS INCORRECT YOU 

MUST RE-ENTER 

830 PRINT "EVERYTHING." 

840 PRINT "" 

850 PRINT "IS ALL DATA CORRECT?" 

860 INPUT HS 

870 IF HS--NO" OR HS-"N" THEN 550 ! REPEATS Spec_data 

880 PRINT — 

890 PRINT "DO YOU WANT A HARD COPY?" 

900 INPUT HS 

910 IF HS-”N" OR HS-”N" THEN 1060 
920 ! 

930 PRINTER IS 706 ! TOGGLES PRINTER TO THERMAL PRINTER 

940 PRINT "" 

950 PRINT "USING ENGLISH UNITS (in.Ibs.psi) THE FOLLOWING INFORMATION” 

960 PRINT "HAS BEEN ENTERED FOR: ";DS 
970 PRINT "" 

980 PRINT "WIDTH -*\-W:" NOTCH LENGTH -”;A 

990 PRINT "Bmax "”:Bmax:" Bmui *";8min 

1000 PRINT "HALF SPAN -":Hs:" POISSON'S RATIO >";Nu. 

1010 PRINT "YIELD STRENGTH -";Ys;" ULTIMATE TS -”;Uts 
1020 PRINT "ELASTIC MODULUS -";E;" FLOW STRESS -";Fs 
1030 PRINT "" 

1040 PRINTER IS 1 ! TOGGLES PRINTER BACK TO CRT 

1050 ! 

1060 PRINT CHRS< 12> ! CLEARS CRT 

1070 ! 

1080 ! CONVERSION TO METRIC UNITS 

1090 Conlen»25.4 ! 1 in - 25.40mm 

1100 Wm”W»Conlen 

1110 Am-A»Conlen 

1120 Bmaxm~Bmax*Conien 

1130 Bminm“Bmin»Conien 

1140 Hsm-Hs**Conlen 

1150 Conpress-6. 89475293 ! tpsi - 6 . 89475293kPa (HP PROGRAMING BOOK) 

1160 Em-£~Conpress 
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1170 Ysm* Ys'Conpress 
1180 Utsm*Uts»Conpress 
1190 Fsm“Fs»Conpress 
1200 1 

1210 PRINT "CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) for ";DS;" YIEL 
DS:" 

1220 PRINT "UIOTH *":Wm:" NOTCH LENGTH »";Am 
1230 PRINT "Bmax *":Bmaxm:" Butin •"iBnuntn 

1240 PRINT "HALF SPAN -”;Hsm:" POISSON'S RATIO -";Nu 

1250 PRINT "YIELD STRENGTH -";Ysm;” UTILMATE TS >":Utsm 

1260 PRINT "ELASTIC MODULUS -" ;En» ;** FLOW STRESS >";Fsm 

1270 PRINT "" 

1280 PRINT "DO YOU WISH TO HAVE A PRINT OUT OF THIS DATA?" 

1290 INPUT G3 

1300 IF GS-"N” CR GS-"NO" THEN 1410 
1310 ! 

1320 PRINTER IS 706 1 TOGGLES TO THERMAL PRINTER 

1330 PRINT "CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) FOR: ":DS 

1340 PRINT "WIDTH >":Wm;" NOTCH LENGTH -" ; An. 

1350 PRINT "Bmax *";8maxm:" Bmin 

1360 PRINT "YIELD STRENGTH «":Ysm;" ULTIMATE TS -";Utsm 

1370 PRINT "ELASTIC MODULUS -";Em;" FLOW STRESS -";Fsm 

1380 PRINT "" 

1390 PRINTER IS 1 ! TOGGLES PRINTER TO CRT 

1400 ! 

1410 PRINT CHRS(t2) ! CLEAR CRT 

1420 ! 

1430 Flag3-0 
1440 A 1 =A 
1450 GOSUB Params 
1460 ! 

1470 RETURN 
1480 ! 

1 490 ! * ***** 

1500 ! 

1510 Neu_a: ! THIS ROUTINE COMPUTES NEW CRACK LENGTH AND NEW PARAMETERS 

1520 ? 

1530 Bef fm*Bmaxm- (Bmaxm-Bminm) '2/Bmaxm 1 B EFFECTIVE FROM ASTM E813-81 

1540 02- 1 » INITIALIZES COUNTER 

1550 PRINT CHRSU2) ! CLEARS CRT 

1560 PRINT "THIS ROUTINE WILL OBTAIN 'HJRRENT CRACK LENGTH- 
1570 PRINT M AND NEW PARAMETERS FOR PRECRACKING." 

1580 PRINT "BEFORE YOU CAN RUN THIS ROUTINE YOU MUST:" 

1590 PRINT " 1. ENTER SPECIMEN DATA" 

1600 PRINT ” 2. ENTER CONVERSION DATA” 

1610 PRINT " 3. SETUP MTS FOR MANUAL LOADING OF THE”; 

1620 PRINT ” SPECIMEN.” 

1630 PRINT " ” 

1640 PRINT " ARE YOU READY?” 

1650 INPUT HS 

1660 IF HS-"N” OR H$-"NO” THEN 2450 * RETURNS TO Ir-KEY SELECTION 

1670 PRINT CHRS<12> ! CLEARS CRT 

1680 PRINT "YOU MAY LOAD AND UNLOAD SPECIMEN AS OFTEN AS YOU LIKE — BUT YOU M 
UST REMAIN 

1690 PRINT "BELOW CURRENT PMAX OF”:Pmaji :"Lbs" 

1700 DIM Lod<2000> .Def K2000) ,Cod<2000) 

1710 DIM Lodm(2000) .Def lm<2000) ,Codm(2000) 

1720 PRINT "WHEN YOU ARE READY TO BEGIN PRESS THE 'ENTER' OR 'CONTINUE' KEY." 
1730 INPUT HS 

1740 PRINT CHRS< 12) ! CLEARS CRT 
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! USED AS COUNTER 
! USED AS INDEX 



K-0 
1*1 
Fla 

Flag2-0 

WHILE F l ag I - 0 ! ESTABLISHES ONE LOAD-UNLOAD DATA TAKING LOOP 

GOSUB Data_acq 
END WHILE 

Flag2*0 ! RESETS FLAG2 

WHILE Flag2*0 

GOSUB Data_acq 
END WHILE 

f 

I THIS SECTION CAN BE USED TO VIEW THE COLLECTED DATA BY REMOVING THE REMA 
RK CHARACTER <!>. 

1880 ! FOR 1-1 TO K 

1890 ! PRINT USING 1 960 ; I .Lodmt I) , Def Im< I ) ,Codm< I) 

1900 1 IMAGE X.3D,4X,6D.D.4X.3D.7D.7X.3D.7D 

1910 ! NEXT I 

1920 ! WAIT 5 

1930 K1-1 

1940 K3-K 

1950 GOSUB Slope 1 SLOPE - I/COD COMPLIANCE 

1960 GOSUB Cracic_length ! USES COD COMPLIANCE 
1970 ! 

1980 AO < D2) -Oracle lm/25 .4 ! CONVERTS TO INCHES 

1990 A1-A0< D2> 

2000 ! 

2010 ! FIND 95% CONFIDENCE INTERVAL FOR CRACK LENGTH 

2020 St i f f m-S t if fm+T9»U4 
2030 GOSUB Oracle length 

2040 Amax-Cracle lm/25 . 4 ! CONVERTS TO INCHES 

2050 Stiffm-Stiffn-2«T9-U4 
2060 GOSUB Crack_length 

2070 Amtn-Cracie lm/25. 4 t CONVERTS TO INCHES 

2080 1 
2090 ! 

2100 PRINT CHR$< 12) * CLEARS CRT 

2110 PRINT "THE FOLLOWING INFORMATION HAS BEEN OBTAINED:" 

2120 PRINT "" 

2130 PRINT "CRACK LENGTH «":A0<D2> :"in FROM MEASUREMENT *";D2 
2140 PRINT "WITH 95Z CONFIDENCE INTERVAL FROM";Amin;" TO" ;Amax 
2150 PRINT 

2160 PRINT "USING COD COMPLIANCE OF:";Compl ;"in/Lb ;Comp Im ; "m«/N" 

2170 PRINT "WITH A CORRELATION OF:";Correl :" FROM" ;K; "DATA POINTS" 

2180 PRINT "" 

2190 PRINT “DO YOU WANT A HARD COPY?" 

2200 INPUT GS 

2210 IF G$-”N” OR GS-"NO" THEN 2340 
2220 ! 

2230 PRINTER IS 706 ! TOGGLES TO THERMAL PRINTER 

2240 PRINT 

2250 PRINT "THE FOLLOWING INFORMATION HAS BEEN OBTAINED:" 

2260 PRINT "CRACK LENGTH »";A0<D2) :"in FROM MEASUREMENT #";D2 
2270 PRINT "WITH 95XC0NFIDENCE INTERVAL FROM" ; An* i n ; " TO";Amax;"in" 

2280 PRINT "USING COD COMPLIANCE OF:" :Compl :"in/Lb -";Complm;"mw/N" 

2290 PRINT "WITH A CORRELATION OF:" :CorreI ;"FROM" ;K;“ DATA POINTS." 

2300 PRINT 

2310 PRINTER IS 1 * TOGGLES PRINTER TO CRT 

2320 ! 



1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1360 

1870 
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2330 PRINT 

2340 PRINT “ DO TOU WANT ANOTHER MEASUREMENT?” 

2350 INPUT FS 

2360 D2-D2+1 * INCREMENTS COUNTER 

2370 IF F$-’*Y" OR FS-“YES“ THEN 1670 
2380 * 

2390 IF A0<D2-1 )>0pta THEN G0SU8 End_precrack ! PRECRACKING IS COMPLETE 
2400 ! 

2410 IF A0< D2- 1 ) > . 4-De la+A THEN Flag3-1 ! SETS FLAG INDICATING IN LAST 607. 0 
F PRECRACKING 
2420 * 

2430 GOSUB Params 
2440 ! 

2450 RETURN 
2460 ! 

2470 ! * «**.**«** 

2480 ! 

2490 Data.acq: ! DATA ACQUISITION ROUTINE 

2500 OUTPUT Dvm:"VR5 AF 1 AL3" 

2510 ? 

2520 OUTPUT Dvm;"AI1 VTf f READS LOAD 

2530 ENTER Dvm;Lod<I) 

2540 Lodm< I)-Lod<I)»$ln»< 1 ) 

2550 ! 

2560 OUTPUT Dvm:'*AI2 VT1'* ! READS COD 

2570 ENTER Dvm:CodU> 

2580 Codrr»(I)-Cod<I)-SIir»<2> 

2590 ! 

2600 OUTPUT Dvn:'*AI3 VTt” ! READS DEFLECTION 

2610 ENTER Dvm:DefKI> 

2620 Def Im(I)-Def KI)-SIn»<3) 

2630 ! 

2640 ! TEST DATA TO BREAK OUT OF DATA ACQ LOOP 

2650 ! 

2660 IF I- I THEN 2720 • ! SKIPS IF FIRST DATA POINT 

2670 Lodif fm”Lodm( I )-Lodm< 1-1 ) 

2680 IF Lodi ffm< -25.0 THEN Flagl-1 * MUST BE UNLOADING 

2690 IF Lodiffm>25.0 THEN Flag2-1 ! MUST BE RELOADING. 

2700 Lodiffm-0 
2710 ! 

2720 I-I+l ! INCREMENTS INDEX 

2730 ! 

2740 K=*K<-1 ♦ INCREMENT COUNTER 

2750 t 
2760 RETURN 
2770 ! 

2780 » ..................w.................*..., 

2790 ! 

2800 Craclc_length: ! CALCULATES CRACK LENGTH USING HUDAK (DTNSRDC PROG) 

2810 ! NOTE: THIS ROUTINE USES COD COMPLIANCE METHOD VICE DEFLECTION COMPLIA 

NCE . 

2820 FOR 1-1 TCI 6 
2820 READ Coef(I) 

2840 NEXT I 

2850 Complm-1 /Stif fm 

2860 C*2»Wm«Bef f m-Em-Complm/ <Hsm» 1000) 

2870 C-1/<SQR<C)+1 ) 

2880 G< 1 )-C 

2890 FOR 1-2 TO 5 

2900 G< I)-G< I- I )"G< 1 ) 
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2910 NEXT I 

2920 F*Coef ( l ) 

2930 FOR 1*1 TO 5 

2940 F*F+Coef( 1*1 >*G(I) 

2950 NEXT I 

2960 Crack ln**F«Wm 

2970 RESTORE 

2980 DATA 0.998265.-3.81662,-1 .80596,32.31038,-44.15665.-52.67376 
2990 ! 

3000 RETURN 

3010 ! 

3020 1 * - 

3030 ! 



3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

2150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3230 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 



Conv_data: ! INPUTS CONVERSION DATA FOR MEASURING DEVICES 

PRINT CHR5< 12) ! CLEARS CRT 

PRINT "INPUT: LOAD RANGE, STRAIN RANGE. STROKE RANGE" 



PRINT 

PRINT " (BY PERCENT ie. 

INPUT Lr.Cr.Sr 
PRINT "" 

PRINT "LOAD RANGE - ":Lr 
PRINT "STRAIN (COD) RANGE 



tO. 20. 50 OR 100)’ 



':Cr 



:Sr 



PRINT "STROKE (DEFLECTION) RANGE ■ 
Lmax-20000 ! MAX LOAD RANGE IN (lbs). 

Cmax* .200 ! MAX COD RANGE IN (tn). 

! MAX STROKE RANGE IN (in). 
1 MAX VOLTAGE IN (volts). 



Smax*5 . 0 
V— 10 



? MAXIMUM OPERATING POINTS ARE: 

0 1 max * Lmax »L r / 1 0 0 

0cmax-Cm*x*Cr/t00*<-1.0) * -1 CHANGES SIGN 

0smax*Snax«Sr/1 00 

; 

I CONVERSION CONSTANTS ARE: (ENGLISH UNITS) 
S l ( 1 )*0 imax/V 
SI (2)*0cmax/V 
Sl(3)*0smax/V 
! CONVERTING TO METRIC: 



! 25 . 4mm * 1 in 
! 4.448221615N • 



1 lb 



Con len»25 . 4 
Conlb-4. 448221615 
Slm(1)*Sl(1 )»Conlb 
Slm(2)*S I (2>»ConIen 
SIm(3)*SI (3)**ConIen 
» 

PRINT ”" 

PRINT 

PRINT "YOUR ENTRIES HAVE PRODUCED THE FOLLOWING 



(HP PROGRAMING BOOK) 



CONVERSION DATA 



3330 PRINT ”•* 

3390 PRINT "FOR LOAD: ”;Slm( l > ;"NEWTON/VOLT (";S1(1 ) :"Lbs/V) 
3400 PRINT "•• 

3410 PRINT "FOR COD: ";SIm(2) :"mm/V (";SI(2> :"in/V)" 

3420 PRINT "" 

3430 PRINT "FOR DEFLECTION: ";Slm(3> ;"mm/V <";SK3) ;”in/V)" 
3440 PRINT "" 

3450 PRINT "WOULD YOU LIKE A HARD COPY?" 

3460 INPUT H$ 

3470 IF H$*"N" OR H5-"N0" THEN 3570 

3480 ! 

3490 PRINTER IS 706 ! TOGGLES PRINTER 
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3500 PRINT — 

3510 PRINT •• CONVERSION DATA" 

3520 PRINT ’’FOR LOAD : " ;Slm< 1 > : "N/V <" ;S1 ( 1 > ; ”Lbs/V >" 

3530 PRINT "FOR COD:" :SIm<2) :"rom/V <":$l<2> :"m/V>" 

3540 PRINT "FOR DEFLECTION:" ;Sln<3> :"mm/V <":S1<3) in/V) 
3550 PRINTER IS 1 ! TOGGLES PRINTER 

3560 ! 

3570 PRINT CHRSU2) 

3580 ! 

3590 RETURN 
3600 1 

3610 * 

3620 ! 



3630 

3640 

3650 

3660 

3670 

3630 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 



Slope: ! LINEAR LEAST SQUARES FIT 

PRINT CHRS< 12) ! CLEAR CRT 

1 INITIALIZE VARIABLES 



Sumx-0 
Sumx2“Q 
Sumy“Q 
Sumy2“0 
Sumxy “0 
FOR I-K1 

f 



1 SUM OF X VALUES 

* SUM OF X SQUARED VALUES 

! SUM OF Y VALUES 

! SUM OF r SQUARED VALUES 

• SUM OF THE PRODUCT 

TO K3 ! KJ-INDEX OF FIRST. K3-INDEX OF LAST DATA POINT USED. 



X”Codm< I ) 

Y“Lodm< I> 

i 

Sumx”Siunx*X 
Sumy “Sumy+Y 
Sumx2*Sumx2+X*X 
Sumy 2 ” S umy 2 *■ Y * Y 
Stunxy“Sujnxy + X»Y 
NEXT I 

K3-K3-P2 ! USED TO CALCULATE NUMBER OF DATA POINTS WITHIN BOUNDS. 

K4-K3-K1+1 ! K4-T0TAL * DATA POINTS USED IN CALC 

Qal “Suax2 - Suax*Sujnx/K4 

Qa2”Sumxy-Sumx“Suniy/K4 

Qa3“Suay2-Suoy*Suffly/K4 



M-Qa2/Qat 



! SLOPE 



IF M<0 THEN 4190 ! IF SLOPE <0, BAD DATA 

FROM DTNSRDC PROGRAM 
Q1 “<K4»Sumx2-Sumx~2)/<K4«<K4-1 ) ) 
Q2“<K4»Sujny2-Siuny*2>/(K4“<K4-1 )) 

Q3-OC4-1 )-<G2-M*rt»Q1 >/<K4-2> 

Q4-SQR(Q3/(Q1»<K4-1 > ) ) 

! 

B”<Sumy-N*SuBx)/K4 * INTERCEPT 

i 

Correl-M»SQR(Q1/Q2> 

* 

Qa4“SQR< A8S< ( Sumy2 ~B“ Sumy _ M“S urn xy ) /(K4-2)) ) 
Sigi»aJ3”Qa4/SQR< A8S<Qa1 ) ) 

Sigmam-S igraab»SQR(Sumx2/K4) 

Stiffm-M ! STIFF - SLOPE 

Complin- 1 /M ! COD COMPLIANCE - 1 /SLOPE 

Comp I -Comp lm» < 4 . 44822 1615) /25 . 4 ! CONVERTS TO 

» 

i STUDENT t TEST 
U4“Q4 



in/LB 
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1 RETURNS 



4'00 T9-!. 36 

4110 IF K4>120 THEN 4170 
4120 T9-2.1 

4130 IF <4<20 THEN 4170 ! RETURNS 

4140 N9-K4-2 

4150 T9-N9/<-.6l7109+N9-. 510239) 

4160 ! 

4170 RETURN 
4130 ! 

4190 PRINT "POOR DATA — SLOPE < 0” 

4200 M-333333.3 

4210 Correi*0 
4220 ! 

4230 RETURN 
4240 ! 

4250 ! * 

4260 ! 

4270 Parans : * THIS ROUTINE WILL DETERMINE THE CURRENT MAX AND MIN LOAD FOR 

FATIGUE PRECRACKING 
4280 ! 

4290 ! NOTE THE USE OF ENGLISH UNITS DUE MTS READOUT BEING IN POUNDS. 

4300 ! INITIALIZE VARIABLES 

4310 Faw*0 

4320 Pmajc-0 

4330 Pmin^O 

4340 Delp*0 

4350 Kmaxp”0 

4360 Kmi.np”Q 

4370 Delk-0 

4380 Delke-0 

4390 ! 

4400 0pta-.6-W ! OPTIMUM CRACK LENGTH 

4410 Dela-Opta-AI ! EXTENSION REQUIRED 

4420 Dela50-Dela-< .5>+A ! CRACK LENGTH AT 50% OF EXTENSION 

4430 * 

4440 Pmax-<4/3>-<Bmin-<W-Al)'2*Fs/<2»Hs)) * FROM ASTM E813 

4450 IF Flag3-1 THEN Pma*-< .4>-Pmax ! MAX LOAD DURING LAST 507. OF PRECR 

ACKING (ASTM E813) 

4468 Pmin”Pmax- .95»Pmax ! FROM ASTM E8I3 

4470 Delp*Pmax-Peun 
4480 ! 

4490 ! DETERMINE Delk BASED ON P USING ASTM E399 

4500 ! 

4510 Faun-3-<Al/W)“.5-<1 .99-<Al/W>«< 1 -Al/U>-<2. 15-3 . 93-A1/U+2 . 7-< A 1/U) "2) ) 

4520 ! 

4530 Faud*2»< 1 +2"A1/W>«< t-Al/W) ~1 .5 
4540 ! 

4550 Fay*Fawn/Fajud 
4560 ! 

4570 ! 

4580 Kmaxp-<Proax«2-Hs/<Bmin-U*1 .5)>-Fau ! FROM ASTM E399 

4590 Kminp”<Pmin«2«Hs/(Biiun*W*1 .5) >»Fajj ! FROM ASTM E399 

4600 Deik*Kn»axp-Kmmp 
4610 ! 

4620 ! DETERMINE Delke BASED ON MODULUS FROM ASTM E813 

4630 IF Flag3-0 THEN Delke-40000 * DURING EARLY EXTENSION 

4640 IF Flag3* 1 THEN Deike-< . 00 1 )-E ! DURING LAST 507. OF PRECRACKING IRE 

F. ASTM E399 4 E813 3 

4650 PRINT 
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PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



"BASED ON CURRENT INFORMATION:" 

»••• 

"CURRENT CRACK LENGTH «";Al :"in" 
"OPTIMUM iC*" ;Gpca:" in" 

"EXTENSION REQUIRED -";De!a 
"CRACK LENGTH AT 507. OF EXTENSION 
'•Pmax “": Pmax : "Lbs Pnin “";Pmin: 
"Kmaxp -’*; Kmaxp : " k'minp -":Kminp 
"DELTA_K »";Delk ;"m ' .5" 

"DELTA.KE -";Delke:"in* .5" 

"FAW -" ;Faw ! SERVES AS A CHECK 



*” ;De la50 ; "in" 
Lbs" 



Pmax 



MUST BE DECREASED’ 
SKIP IF OKAY 



*IF Delb > Delke. 

IF Delk<Delke THEN 4900 
! COMPUTE NEW Pmax 
Deipl *<Delke/Fau/<2“Hs)»Bmin»W~1 .5) 
Pmax 1 “Delpl / .95 
Pminl ’Pmax 1 - .95 “Pmax 1 
Kmax I • < Pmax I »2*Hs/ < Bm i n*W 1 . 5 )) »Fau 
PRINT "THE NEW Pmax BASED ON Delke IS: 
•NEW Pm in -";Pmml :"Lbs" 

•NEW DELTA_P -";Delp1 ;"Lbs" 

•NEW Kmax -";Kmaxl 



DO YOU WANT A HARD COPY?* 



PRINT 
PRINT 
PRINT 
» 

PRINT " 

INPUT HS 

IF HS-”N" OR HS--NO" THEN 5080 
PRINTER IS 70S ! TOGGLES PRINTER 

PRINT "BASED ON CURRENT INFORMATION FOR 
PRINT "" 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 



1 ;Pmax 1 ; "Lbs' 



' ; DS ; — : 



4860 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 

5180 ! 

5190 A_table : 

S CRACK LENGTHS 
5200 ! 

5210 PRINT "BEFORE YOU CAN RUN THIS SUBROUTINE YOU MUST ENTER:’ 
5220 PRINT" 1. SPECIMEN DATA” 

5230 PRINT " 2. CONVERSION DATA" 



CURRENT CRACK LENGTH -":At ;"in" 

OPTIMUM CRACK LENGTH ;Opta:"in" 

EXTENSION REQUIRED -":Dela:"in" 

CRACK LENGTH AT 502 OF EXTENSION «";Dela50 ;"ir 
Pmax »’*;Pmax;’*Lbs Pm in -":Pmin;’*Lbs" 

Delp “”;Delp :’*Lbs Kmax *":Kmaxp 
DELTAJC *’*;Deik :’*in* .5 DELTA_KE -";Delke;’*in - 
IF Deik<Delke THEN 5090 
PRINT 

NEW LOADS BASED ON DELTA.KE ARE:" 

Pmax Pmax 1 : "Lbs Pm in •" ;Pm ini ; "Lbs" 
DELTA_P “”;DeI P 1 ;**Lbs Kmax -";Kmax1 



.5" 



PRINT 
PRINT 
PRINT 
PRINT 
PRINTER 
PRINT 
PRINT 
PRINT 
INPUT 

I 

RETURN 

t 



IS 1 



! TOGGLES PRINTER TO CRT 



"MAKE THE NECESSARY ADJUSTMENTS BEFORE CONTINUING.’ 
"WHEN COMPLETED PRESS 'CONTINUE'." 

CS 



! PRODUCES TABULATED VALUES FOR LOAD vs COD VALUES AT VARIOU 
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5240 PRINT 

5250 PRINT "ARE YOU READY - )" 

5260 INPUT HS 

5270 IF HS-"N" OR H$-"NO" THEN 5550 ! RETURNS TO Ic-KEY SELCETION 

5280 ! 

5290 BefF-Bmax-<Bmax-Bfiun)*2/Bmax ! IAW ASTM E8I3-31 
5300 PRINT "BEFF -”;Beff 
5310 ! 

5320 ! USING CRACK LENGTH EQUATION FROM DTNSRDC PROGRAM CALC CRACK. LENGTH FOR 

VARIOUS VALUES OF U 

5330 A00-. 998265 

5340 A1— 3.81662 

5350 A2— 1.80596 

5360 A3-32.31038 

5370 A4— 44.15665 

5380 A5--52. 67876 

5390 FOR U» .07 TO .155 STEP .005 

5400 Ar-A00«-A1 »U+A2-U*2+A3-U'3+A4*U*4+A5-U*5 

5410 CH-Ar-H 

5420 PRINTER IS 706 ! TOGGLES TO THERMAL PRINTER 

5430 PRINT "U -”;U 

5440 PRINT "a/W -":Ar:"in CRACK LENGTH -" :Ck 1 : "in" 

5450 PRINT "LOAD (lbs) COD (in) COD (Volts)" 

5460 FOR Lodi -200 TO 1800 STEP 100 

5470 C-Bef f •H*E*2/Hs 

5480 Cod 1 -Lodi »< < 1 -U)/U) ~2/C 

5490 Codv-Cod1/Sl<2) 

5500 PRINT Lodi , Cod l ,Codv 

5510 NEXT Lodi 

5520 NEXT U 

5530 PRINTER IS 1 ! RETURNS TO USING CRT AS PRINTER 

5540 ! 

5550 RETURN 
5560 ! 

5570 * — 

5580 ! 

5590 End_precracK: ! END OF PRECRACKING 

5600 ! 

5610 PRINT “ PRECRACKING IS COMPLETED!" 

5620 PRINT "" 

5630 PRINT "FINAL CRACK LENGTH -”;A0<D2-t> 

5640 PRINT "CORRELATION OF DATA -";Correl 

5650 PRINT "" 

5660 PRINT " REMOVE SPECIMEN FROM MTS." 

5670 ! 

5680 END 
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APPENDIX E 



J INTGL LISTING 



10 

20 

30 

40 

50 

60 

70 

80 

90 

too 
I to 
120 
130 
140 
150 
160 
170 
180 
130 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
530 
590 



! J-INTEGPAL FRACTURE TOUGHNESS PROGRAM 

* FOR THREE POINT BEND SPECIMEN. THIS 

* PROGRAM INCORPORATES THE USE OF THE 
! FOLLOWING EQUIPTMENT: 

! t . HP-9826 COMPUTER 
! 2. HP -3 49 7 A OATA ACQUISITION SYSTEM 

! 2. HP-3437A SYSTEM VOLTMETER 

! 4. HP-2671 G THERMAL PRINTER 

! 5. HP-7225B PLOTTER 

! THIS PROGRAM WAS WRITTEN BY: 

? LCDR H . K . TRITCHLER. USN 
! VERSION 1.0 DATE 11/33 



COM /J_intgi/ Fl 
» 

I DIN Lod'8000) .Def K8000) .Cod(3000> 

! NOTE: Lod. Def I AND Cod ARE ACTUALLY DIMENSIONED IN Imtiai_a ROUTINE 
? Lod< I ) - LOAD 
! Def 1(1) - DEFLECTION 

! Coa<I> » CRACK OPENING DISPLACEMENT (COD) 



DIN AO (50) .Amin0(50> .Ar»ax0(50) .Como 1 0(50) ,Correi0(50> 
DIM Area9( 100) .AK100) .Coef (6) .G<6) 



DIM J9( 100) . Delta_a( 100) ,Mdef 1(100) 
DIM DataS(1 :4)C10] 

DUMP DEVICE IS 706 
Dvm*709 
CLEAR Dvm 
PRINT CHRS( 12) 



! THERMAL PRINTER ADDRESS 
! DATA ACQUISITION ADDRESS 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



! CLEARS CRT 

- J- INTEGRAL TEST PROGRAM" 
VERSION 1.0" 

WRITTEN BY LCDR W.K. TRITCHLER" 
PLEASE SELECT A 1c -KEY" 



Begin: ! 

ON KEY 0 LABEL "SPEC_DATA" GOSUB Spec_data 
ON KEY 1 LABEL "CONV_DATA" GOSUB Conv_data 
ON KEY 2 LABEL "DAT_FILE“ GOSUB Data_file 
ON KEY 3 LABEL "PLOTTER" GOSUB Setup plotter 
ON KEY 4 LABEL "GCLEAR" GOSUB Gclr 
ON KEY 5 LABEL "aO" GOSUB Initial_a 
ON KEY 7 LABEL "START" GOSUB Main _prog 
!<-) ON KEY 3 LABEL "R-CURVE" GOSUB R_curve 
! (SEE NOTES ON («> IN R_curve ROUTINE 
ON KEY 9.4 CALL Stop_test 



Spin: ! WAITS FOR KEY INPUT 

GOTO Spin 
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600 Gel r : ! GRAPHICS SCREEN CLEAR 

610 ! THIS ROUTINE ALLOWS SINGLE KEY STROKE TO REMOVE A GRAPHICS DISPLAY C RCM 

CRT 

620 GCLEAR 

630 RETURN ! RETURNS TO It -KEY SELECTION 

640 ! 

650 !** — 

650 ! 

670 Spec.data: ! INPUTS REQUIRED SPECIMEN DATA 

660 PRINT CHRSM2) ! CLEARS CRT 

630 PRINT "INPUT SPECIMEN IDENTIFICATION (ID):" 

700 INPUT DS 

7?0 PRINT ”” 

720 PRINT "J-INTEGRAL TEST for ";DS 

730 PRINT "" 

740 PRINT "INPUT SPECIMEN DATA USING ENGLISH UNITS (LBS. INCHES, etc.)" 

750 PRINT "" 

760 PRINT "WHAT IS THE SPECIMEN'S WIDTH. NOTCH LENGTH, Bmax . 8mm (THICKNESS) 
?" 



770 INPUT W, A, Bmax. Bmm 

7B0 PRINT "" 

790 PRINT "WHAT IS THE SPECIMEN'S HALF SPAN AND POISSON'S RATIO?” 

800 INPUT Hs.Nu 

8 * 0 PRINT "" 

820 PRINT "WHAT IS THE SPECIMEN'S ELASTIC MODULUS. YIELD STRENGTH AND ULT 

IMATE STRENGTH?" 



830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
MUST 
950 
960 
970 
980 
990 
1000 
1 0 1 0 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1 120 
1130 
1140 
1150 



INPUT E.Ys.Uts 

PRINT CHRS( 12) ! CLEARS CRT 

Fs-< Ys+Uts)/2 

PRINT "USING ENGLISH UNITS ( in . lbs . psi )" 

•THE FOLLOWING HAS BEEN ENTERED for: 
'WIDTH - " ; W : “ NOTCH LENGTH -”;A 

""Bmax »”;Bmax:” Bmm *”;8mm 

•HALF SPAN -”:Hs;" POISSON'S RATIO 



■;DS 



•YIELD STRENGTH 
•ELASTIC MODULUS 



•;Ys:' 
■ * ; E : * 



ULTIMATE TS 



' ;Nu 
’ :Uts 



FLOW STRESS -";Fs 
THIS PROGRAM ALLOWS NO EASY RECOVERY, IF ANY DATA IS INCORRECT YOU 



! REPEATS Spec_data 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
RE-ENTER 

PRINT "EVERYTHING." 

PRINT "" 

PRINT "IS ALL DATA CORRECT?" 

INPUT HS 

IF HS-"NO" OR HS-"N" THEN 670 
PRINT ”" 

PRINT "DO YOU WANT A HARD COPY?" 

INPUT HS 

IF HS-"N" OR HS-”N" THEN 1180 
» 

PRINTER IS 706 ! TOGGLES PRINTER TO THERMAL PRINTER 

PRINT "" 

PRINT "USING ENGLISH UNITS (m.lbs.psi) THE FOLLOWING INFORMATION' 

PRINT "HAS BEEN ENTERED FOR: ";DS 

PRINT 

PRINT "WIDTH -" :W:" NOTCH LENGTH «";A 

PRINT "Bmax *":Bmax;" Bmm »";Bmm 

PRINT "HALF SPAN POISSON'S RATIO »";Nu 

PRINT "YIELD STRENGTH -";Ys:" ULTIMATE TS -";Uts 
PRINT "ELASTIC MODULUS -”;E:” FLOW STRESS -":Fs 
PRINT "" 
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1 1 60 PRINTER IS 1 ! TOGGLES PRINTER BACK TO CRT 

1170 I 

t i SO PRINT CHR3(12) ! CLEARS CRT 

1190 ! 

1200 * CONVERSION TO METRIC UNITS 

1210 Conlen*25.4 ! 1 m * 25.40mm 

1220 W-U-Conlen 

1230 A»A»Conlen 

1240 Bmax*Bmajc»Con i en 

1250 Bmin-Bmin»Conien 

12S0 Hs*Hs*Conlen 

1270 Conpress-6. 39475293 * Ipsi - 6 .39475293kPa (HP PROGRAMING BOOK) 

1280 E-E»Conpress 
1290 Ys*Ys*Conpress 
1300 Uts’Uts«Conpres3 
1310 Fs*Fs»Conpress 
1320 ! 

1330 PRINT "CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa/ for ";DS:" YIEL 
DS:" 

1340 PRINT "WIDTH -";W:" NOTCH LENGTH *":A 

1350 PRINT "Bmax -”:Bmax:" Bmtn 

1360 PRINT "HALF SPAN POISSON'S RATIO -";Nu 

1370 PRINT "YIELD STRENGTH «";Ys:" UTILMATE TS • M ;Uts 

1380 PRINT "ELASTIC MODULUS -";£;" FLOW STRESS -";Fs 

1390 PRINT "" 

1400 PRINT "DO YOU WISH TO HAVE A PRINT OUT OF THIS DATA?” 

1410 INPUT GS 

1420 IF GS«"N” OR G3>"N0" THEN 1540 
1430 * 

1440 PRINTER IS 706 ! TOGGLES TO THERMAL PRINTER 

1450 PRINT "CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) FOR: ";DS 

T460 PRINT "WIDTH -”:W:" NOTCH LENGTH -";A 

1470 PRINT "Bmax -":Bmax:" Bmtn »":Bmtn 

1480 PRINT “HALF SPAN -":Hs:” POISSON'S RATIO -":Nu 

1490 PRINT "YIELD STRENGTH -”;Ys;" ULTIMATE TS *”;Uts 

1500 PRINT "ELASTIC MODULUS FLOW STRESS -";Fs 

1510 PRINT ”” 

1520 PRINTER IS 1 ! TOGGLES PRINTER TO CRT 

1530 ! 

1540 PRINT CHR$< 12) ! CLEAR CRT 

1550 ! 

1560 RETURN ! RETURNS TO It -KEY SELECTION 

1570 ! 

1580 ! 

1530 ! 

1600 Setup_plotter : ! SET UP INTERACTIVE PLOTTER 
1610 ! YOU CAN CHOOSE BETWEEN THREE DIFFERENT PLOTS 

1620 ? <P vs COD. P vs DEFLECTION. OR J vs DELTA_A) 

1630 PRINT CHRS( 1 2) * CLEARS CRT 

1640 DIM TST801 , XSI801 , YSI80] 

1650 GCLEAR 
t660 GINIT 
1670 GRAPHICS ON 

1680 PLOTTER IS 3. "INTERNAL” ! PLOTS ON CRT FOR VIEWING BEFORE SENDING TO 

HP-7225B PLOTTER 
1690 PRINT "TYPE OF PLOT:” 

1700 PRINT ” 1. LOAD vs COD" 

1710 PRINT ” 2. LOAD vs DEFLECTION" 

1720 PRINT " 3. J vs CRACK EXTENSION- 

1730 PRINT "ENTER THE * OF THE PLOT YOU DESIRE" 
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1740 

1750 

1760 

1770 

1780 

1790 

1800 

1310 

1820 

1830 

1840 

1350 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 



INPUT C2 
PRINT ,,M 
PRINT "" 

IF C2-1 THEN 1810 
IF C2*2 THEN 1860 
IF C2 = 3 THEN 1910 
PRINT "" 

TS=”PL0T of LOAD vs COD for "ADS 
PRINT TS 
XS»"C0D <mm)" 

YS=”L0AD < N ) " 

GOTO 1950 

TS--PL0T of LOAD vs DEFLECTION for "ADS 
PRINT TS 

XS-"DEFLECTI0N (mm)” 

Y$="L0AD <N>" 

GOTO 1950 

TS--PL0T of J vs CRACK EXTENSION for "ADS 
PRINT TS 

XS«"DELTA_a (mm)" 

YS»”J (Ic J/m'2)“ 

PRINT "" 

PRINT "USING METRIC UNITS (mm.N.k J/m '2. etc) :" 
PRINT "ENTER HERE THE Xmin ,Xmax . Ynun . Ymax" 

INPUT XKC2) ,X2(C2>.Y1(C2> .Y2<C2) 

X3=X2(C2>-X1 (C2) 

Y3=Y2(C2)-Y1 <C2> 

VIEWPORT 13.100-RATI0. 10, 90 
WINDOW XI (C2) ,X2(C2> . Y1 (C2> , Y2(C2> 

PRINT "" 

PRINT "INPUT DELTA_X AND DELTA_Y LABEL INTERVAL" 
PRINT "FOR THE TWO AXES:" 

INPUT X4.Y4 
PRINT "" 

PRINT "DO YOU WISH- TO HAVE A GRID OVERLAY?” 

INPUT GS 

IF G$»"N0" OR G$*"N" THEN 2140 
PRINT "" 

PRINT "INPUT GRID LINES PER DELTA LABEL:" 

INPUT B0 

PRINT "PRESS ’ENTER' TO PLOT" 

INPUT H$ 

PRINT CHRS(12> ! CLEARS CRT < 

IF GS»”N0" OR GS»"N” THEN 2310 
FOR I-X4/B0+X1 (C2> TO X2(C2) STEP X4/B0 
MOVE I.Y1(C2> 

DRAW I , Y2(C2) 

NEXT I 

FOR I-Y4/B0+Y1 <C2) TO Y2(C2) STEP Y4/B0 
MOVE XI <C2) , I 
DRAW X2(C2> , I 
NEXT I 
GOTO 2320 

IF G$-"Y" OR G$»”YES" THEN 2320 
MOVE X2<C2> . Y1 <C2> 

DRAW X2(C2) ,Y2(C2> 

DRAW X1<C2).Y2<C2> 

CLIP ON 

AXES X4 , Y4 ,X 1 (C2> . Y1 <C2> 

CLIP OFF 
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2340 

2350 

22G0 

2370 

2330 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2630 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2390 

2900 

2910 

2920 

2930 



VIEWPORT 13. !00«RAri0.5.90 
WINDOW XI <C2) .X2<C2> ,Y1 <C2) , Y2<C2> 
LORG 4 

CSIZE 3.. 6 ! LABELS X DIGITS 

FOR X*X1<C2> TO X2(C2) STEP X4 
CLIP OFF 
MOVE X.YKC2) 

CLIP OFF 

LABEL USING "K”;X 
NEXT X 

VIEWPORT 4. 100»RATIO, 10,90 
WINDOW XI <C2> , X2 < C2 ) . Y 1 ( C2 > . Y2 ( C2 > 
LORG 2 * LABELS Y DIGITS 

FOR Y = Y 1 <C2) TO Y2(C2> STEP Y4 
CLIP OFF 
MOVE X1<C2>.Y 
CLIP OFF 

LABEL USING ”K” ; Y 
NEXT Y 

VIEWPORT O.RATIO*100.0. 100 
WINDOW XI <C2) ,X2<C2) . Y1 <C2) ,Y2(C2> 
CSIZE 5 

LORG 6 ! PLOTS TITLE 

MOVE X2<C2)/2.Y2(C2> 

LABEL TS 



DEG 

LDIR 90 
CSIZE 4.. 7 

MOVE XI <C2) ,Y2(C2)/2 
LABEL YS ! PLOTS Y-AXIS LABEL 

LORG 4 
LDIR 0 

MOVE X2(C2)/2,Y1(C2) 

LABEL XS ! PLOTS X-AXIS LABEL 



PRINT MM 
PRINT " M 
PRINT •"* 
PRINT 
PRINT ” 
PRINT “ M 
PRINT ” 
PRINT " 
PRINT " 
INPUT FS 
IF FS*”N” 



IS THE PLOT (AXES AND LABELS) OK?" 

NOTE: IF YOU ANSWER YES OR Y M 

PLOT WILL BE PRODUCED” 

ON EXTERNAL PLOTTER." 

OR F$-”NO” THEN GOSUB Setup_plotter ! REDO PLOTTER SETUP 



GCLEAR ! CLEARS PLOT FROM CRT 



GOSUB Ext_plotter ! SENDS PLOT TO PLOTTER 

PENUP 
MOVE 0.0 

f 

RETURN ! RETURNS TO Ic-KEY SELECTION OR R.curve ROUTINE 

» 



Ext_plotter : ! DUMPS PLOTTER SETUP TO EXTERNAL HP-7225B PLOTTER 
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ZPmu 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3' 90 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3S10 

3520 

3530 



"'LjT TER iS 705 . “HPCL" 

VIEWPORT 13. 100-RATIG. 10.90 
WINDOW XI <C2> .X2<C2> ,Y1 <C2> . 72 (CO 
IF GS-’W OR GS-”N" THEN 3030 
FOR I-X4/B0+XKC2) TO X2<C2> STEP X4/B0 
MOVE I.TKC2) 

DRAW I.Y2<C2> 

NEXT I 

FOR I* Y4/B0+-Y 1 <C2) TO Y2<C2> STEP Y4/B0 
MOVE X 1 < C2 ) , I 
DRAW X2(C2) , I 
NEXT I 
GOTO 3120 

IF GS--Y" OR G3--YES” THEN 3120 
MOVE X2<C2),YKC2> 

DRAW X2(C2) . Y2(C2) 

DRAW XKC2) ,Y2<C2> 

CLIP ON 

AXES X4.Y4.X1 (C2> , Y1 <C2> 

CLIP OFF 

VIEWPORT t3. 100-RATI0.5.90 
WINDOW XI <C2) .X2<C2> , Y1 <C2) . Y2(C2) 

LORG 4 

CSIZE 3.. 6 ! SIZE OF DIGITS 

FOR X-XKC2) TO X2<C2> STEP X4 
CLIP OFF 
MOVE X.YKC2) 

CLIP OFF 

LABEL USING "K”;X J LABELS X-AXIS 
NEXT X 

VIEWPORT 4. 100-RATIO. 10.90 
WINDOW XI <C2> ,X2<C2> ,Y1 (C2) , Y2(C2) 

LORG 2 

FOR Y»Y1 <C2) TO Y2<C2) STEP Y4 
CLIP OFF 
MOVE XI <C2) , Y ' 

CLIP OFF 

LABEL USING ”K”;Y ! LABELS V-AXIS 

NEXT Y 

VIEWPORT O.IOO-RATIO, 0.100 
WINDOW X1<C2) ,X2(C2) ,Y1 (C2) .Y2(C2) 

CSIZE 5 ! TITLE SIZE 

LORG 6 

MOVE X2(C2)/2.Y2<C2) 

LABEL TS ! PRINTS TITLE 

DEG 

LDIR 90 
CSIZE 4,. 7 

MOVE XI <C2> ,Y2(C2)/2 

LABEL Y$ ! PRINTS Y-AXIS TITLE 

LORG 4 
LDIR 0 

MOVE X2(C2)/2.Y1(C2) 

LABEL XS * PRINTS X-AXIS TITLE 

PENUP 
MOVE 0,0 

t 

RETURN' ! RETURNS TO Setup_pLotter 
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;?.40 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 



Initial 



ESTIMATES INITIAL CRACK LENGTH 



Ber T-8max - < 8max-8m i n ) 2/8max ! USES 3<EFrECTIVE> FROM ASTM E 3 T 3-31 

02- I ! INITIALIZES COUNTER 

PRINT CHRSM2) * CLEARS CRT 

PRINT "BEFORE YOU CAN RUN THIS ROUTINE YOU MUST:" 

PRINT " 1. ENTER SPECIMEN DATA" 

ENTER CONVERSION DATA" . 

SET UP PLOTTER" 

CREATE DATA FILE" 



2 . 

3. 

4. 



PRINT " 

PRINT " 

PRINT " 

PRINT "" 

PRINT "ARE YOU READY?" 

INPUT HS 

IF HS-"N" OR HS»"NO" THEN 5070 ! RETURNS TO KEY SELECTION SCREEN 

PRINT CHRS( 12) ! CLEARS CRT 

PRINT "THIS ROUTINE HILL OBTAIN THE INITIAL CRACK LENGTH AFTER FATIGUE PRE 
-CRACKING." 

3700 PRINT "" 

37 1 0 PRINT "YOU MAY LOAD AND UNLOAD THE SPECIMEN AS OFTEN AS 
OU MUST REMAIN "; 

3720 PRINT "ON THE LINEAR PORTION OF THE CURVE!" 

PRINT "" 

THE MTS SHOULD BE IN 'STROKE' CONTROL" 

AT A VERY SLOW INVERTED RAMP FUNCTION" 

RATE OR YOU MAY DO THIS WITH MANUAL LOAD" 

CONTROL." 



YOU LIKE — BUT Y 



3730 

3740 

3750 

3760 

3770 

3730 

3790 

3800 

3810 

3320 

3830 

3840 

3850 

3360 

3870 

3880 

3P.30 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 



•NOTE : 



‘WOULD YOU LIKE AN ESTIMATE OF Pmax <Lbs>?’ 



OR ES-"NO" THEN 4060 ! SKIPS Pmax CALC 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

f 

PRINT 
INPUT ES 
IF ES-"N" 

f 

PRINT "YOU MUST PROVIDE AN ESTIMATE OF CURRENT CRACK LENGTH- 

PRINT "INPUT YOUR ESTIMATE IN INCHES." 

INPUT Aest 

! CONVERT ALL DIMENSIONS BACK TO ENGLISH UNITS. 

! aO IS SUFFIX ON ALL VARIABLES FOR THIS CALC. 

Bef f a0-8ef F/Conlen 
WaO-W/Conien 

Blig-<WaO-Aest> ! ESTIMATED LENGTH OF UNBROKEN LIGAMENT IN INCHES 

FsaCl-Fs/Conpress $ 

HsaO-Hs/Conlen 

P LaO- .533 “Bef f a0»8 1 ig*2-FsaG/ <2-HsaO) 

•AN ESTIMATE OF Pmax IS":PlaO 
*Pmax ESTIMATES THE POINT AT WHICH THE CURVE" 

•BECOMES NONLINEAR, IAW ASTM E813-81" 



PRINT 
PRINT 
PRINT 

i 

PRINTER IS 706 
PRINT "" 

PRINT "GIVEN a- 
PRINT " Pmax* 

PRINT "" 

PRINTER IS 1 
! 

DIM Lod<8000) ,Def 1(8000) ,Cod(8000) 

f 

! ASSIGN DEFL CORRECTION COEF USED IN Data_acq AND Take_data ROUTINES 
‘ THIS USES A 4th ORDER POLYNOMIAL FIT OF RECORDED DATA ON THE MTS SYSTEM' 



! TOGGLES PRINTER 

Aest: "inches THEN:" 

’ ;PlaO;"Lbs FROM ASTM E813-8V 



S COMPLIANCE 
4100 ! 
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4 1 i u lo c r t I > -5 . 798E-5 
4120 Corr<2>-4 . 7419E-6 

4130 Corr<3)--2.803E-9 



4140 Cor r < 4 ) » 1 . 47 1 E- 1 2 
4150 Corr<5)— 2.774IE-I6 
4160 * 

4170 PRINT •••• 

4180 PRINT "WHEN YOU ARE READY TO BEGIN PRESS THE ENTER' KEY." 

4190 INPUT HS 

4200 PRINT CHRS(12> f CLEARS CRT 

4210 K *0 * USED AS A COUNTER 

4220 I- I » USED AS INDEX 

4230 Flag 1-0 ! NOTE: FLAG1 AND FLAG2 ARE RESET 

4240 FIag2-0 * IN Data_aco ROUTINE 

4250 WHILE Flagl-0 ! ESTABLISHES ONE LOAD-UNLOAD DATA TAKING LOOP 
4260 GOSUB Data_acq 

4270 END WHILE 

4280 FIag2*0 ! RESETS FLAG2 

4290 WHILE Flag2-0 

4300 GOSUB Data_acq 

4310 END WHILE 

4320 PENUP 

4330 ! THIS SECTION CAN BE USED TO VIEW THE COLLECTED DATA BY REMOVING THE REM 

ARK CHARACTER <!>. 

4340 ! FOR 1-1 TO K 

4350 ! PRINT USING 3780 : 1 . Lod< I ) . Def 1 < I) .Cod< I) 

4360 ! IMAGE X ,3D,4X ,6D .0 ,4X .3D . 7D , 7X , 3D . 7D 

4370 ! NEXT I 

4380 ! WAIT 5 

4390 Kl-t 

4400 K3-K 

4410 GOSUB Slope ! SLOPE > I /COD COMPLIANCE 

4420 Cor re 10 (D2) -Cor re I ! SAVES CURRENT CORRELATION COEFFICIENT IN AN ARRAY 

4430 ! 



4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

NT 

4630 

4640 

ft 



Comp 1 0 < D2 > -Comp I 
» 

GOSUB Crack_length 

; 

A0(D2)-Crck 1 

; 

! FIND 957. CONFIDENCE INTERVAL FOR aO 
St if f -Stiff+T9»U4 



! SAVES CURRENT COMPL INACE VALUE IN AN ARRAY 



! USES COD COMPLIANCE 

.e, 

! SAVE THE LAST aO MEASUREMENT OBTAINED IN AN ARRAY 



GOSUB Crack_length 
AmaxO<D2) -Crck 1 
Stiff-Stiff-2«T9*U4 



GOSUB Crack _length 
AminO<D2)-Crck 1 



PRINT CHR$( 12) ! CLEARS CRT 

PRINT "THE FOLLOWING INFORMATION HAS BEEN OBTAINED:" 

PRINT "" 

PRINT "INITIAL CRACK LENGTH <a0) -":A0<D2) ;"mm FROM MEASUREME 

;D2 

PRINT "" 

PRINT "WITH 95% CONDIFENCE INTERVAL ON a FROM"; Ami nO<D2> TO";AmaxO<D2> ;" 



4650 PRINT "" 

4660 PRINT "USING COD COMPLIANCE OF : " ;ComplO<D2> ;"mm/N" 

4670 PRINT "WITH A CORRELATION OF :" ;CorrelO<D2) : "FROM” ;K ; "DATA POINTS." 
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4630 PRINT 

4890 PRINT "DO YOU WANT A HARD COPY?" 

4700 INPUT HS 

4710 IF H$-"N" OR H$-”NO" THEN 4830 
4720 • 

4730 PRINTER IS 706 ! TOGGLES PRINTER 

4740 PRINT "" 

4750 PRINT "THE FOLLOWING INFORMATION HAS BEEN OBTAINED:" 

4760 PRINT "INITIAL CRACK LENGTH *":A0(D2) :"mm FROM MEASUREMENT *";D2 

4770 PRINT "WITH 957. CONFIDENCE INTERVAL FR0M”:Aniin0<D2> TO" ;AmaxO< D2) : "mm" 

4780 PRINT "USING COD COMPLIANCE OF :Comp 10< D2) : "rom/N" 

4790 PRINT "WITH A CORRELATION OF : " ;Corr»iO< D2) : "FROM" ;K ;”DATA POINTS." 

4800 PRINT "" 

4810 PRINTER IS 1 ! TOGGLES PRINTER 

4320 ! 

4830 D2-D2+1 ! INCREMENT COUNTER 

4840 PRINT "" 

4850 PRINT "DO YOU WANT ANOTHER aO MEASUREMENT?" 

4860 INPUT G3 

4870 IF G3«"Y" OR GS-"YES" THEN 4060 
4830 ! 

4890 PRINT "INPUT THE NUMBER <*> OF THE aO MEASUREMENT YOU WISH TO USE AND S 
TORE." 

4900 PRINT " NOTE: ONLY THIS SET OF DATA ARE" 

4910 PRINT " STORED ON DISK." 

4920 INPUT D2 
4930 f 

4940 At »AOf D2) ! ASSIGNS At TO DESIRED A0<D2> VALUE 

4950 ! NOTE: At IS USED IN Main_prog AS INITIAL CRACK LENGTH VALUE 

4960 ! 

4970 PRINT "” 

4980 PRINT "YOU HAVE SELECTED MEASUREMENT *";D2 
49S0 PRINT " ARE YOU SURE?" 

5000 INPUT PS 

5010 ! NOTE: ONLY THIS aO MEASUREMENT IS STORED ON DISK 

5020 ! 

5030 IF PS-"N" OR PS-"NO" THEN 4890 ! ALLOWS CORRECTION 

5040 K2-0 ! K2 IS UNLOADING COUNTER THIS WILL ENSURE CORRECT DAT 

A IS STORED 

5050 GOSvJB Stor_data ! < WRITES ON DISK 

5060 ! 

5070 PRINT CHR$< 12) 

5080 ! 

5090 ! 

5100 RETURN ! RETURNS TO Ic-KEY SELECTION 

5110 ! 

5 f 20 t »**«*««««#*»#*«•«***«««»»•«•**«**«««••••«»* 

5130 ! 

5140 Data_acq : ! DATA ACQUISITION ROUTINE 

5150 OUTPUT Dvm;”VR5 AF1 AL3" 

5160 ! 

5170 OUTPUT Dvm : "AI 1 VT1" ! READS LOAD 

5180 ENTER Dvm:Lod< I> 

5190 Lod< I) -Lod< I)»SI < 1 ) 

5200 ! 

5210 OUTPUT Dvm;”AI2 VT1" ! READS COD 

5220 ENTER Dvm:Cod<I) 

5230 Cod < I ) - < Cod <I))«SI<2) 

5240 ! NOTE: COD IS CALIBRATED ONLY FORM 0.475 TO 0.625 INCHES (MAX COD IS 0.1 

5 INCHES) 
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! READS DEFLECTION 



5250 ! 

5260 OUTPUT Dvn:"AI3 VTJ" 

5270 ENTER Dvm:Defl<I> 

5230 Defl<I>-0efl<I>-Sl<3> 

5290 t 

5300 GOSUB Correct_def l ! DETERMINE DEFL CORRECTION FACTOR 

5310 ' 

5320 Def l < I) -Def l < I) - To tcor 

5330 * 

5340 ! 

5350 ! PLOT THE DATA POINT 

5360 VIEWPORT 1 3 . 1 0O-RATIO . 1 0 .90 

5370 WINDOW X! <C2) .X2<C2> . Y1 <C2) . Y2<C2> 

5380 IF C2-2 THEN 5410 

5390 PLOT Cod< I) ,Lod< I ) 

5400 GOTO 5420 

5410 PLOT Def 1 < I> ,Lod< I > 

5420 ! 

5430 ! TAKE DATA FOR ONE COMPLETE LOADING AND UNLOADING 

5440 IF 1-1 THEN 5520 ! SKIPS TEST ON THE FIRST DATA POINT 

5450 * 

5460 ! TEST DATA 

5470 Lod i f f *Lod< I) -Lod< 1-1 > 

5480 IF Lodif f <-25.0 THEN Flagl-1 * MUST BE UNLOADING 

54S0 IF Lodi ff >25.0 THEN Flag2-1 ! MUST BE LOADING 

5500 Lodiff-0. 

5510 • 

5520 I- 1+1 * INCREMENTS INDEX 

5530 ! 

5540 K-K + t ! INCREMENT COUNTER 

5550 ! 

5560 RETURN ! RETURNS TO Imtial_a ROUTINE 
5570 ! 

5580 ! ******* 

5590 Correct_def I : ! THIS ROUTINE CALCULATES THE DEFLECTION CORRECTION FACT 

OR 

5600 ! 

5610 ! CORRECTION FACTOR IS REQUIRED TO ACCOUNT FOR SYSTEM COMPLIANCE 

5620 ! A 4TH ORDER POLYNOMIAL FIT WAS PREVIOUSLY DETERMINED. 

5630 ! t 

5640 ! NOTE: - — * 

5650 ! THE 5 CORRECTION COEFFICIENTS CORR<I> WERE ASSIGNED IN Initial.a ROUTIN 

E TO SPEED UP PROG. 

5660 ! ........ ............ 

5670 * 

5680 Lodp-Lod< I)/Conib ! CONVERTS CURRENT LOAD TO Lbs FOR CALC 
5690 Def lc< 1 )-Lodp 
5700 FOR J-2 TO 4 

5710 Deflc(J)-Defic<J-1> »Def lc< 1 > 

5720 NEXT J 

5730 Totcor-CorrU) ! PREVIOUSLY ASSIGNED AS NOTED. 

5740 FOR J-1 TO 4 

5750 Totcor-Totcor+Corr< J + 1 ) »Def Lc( J) 

5760 NEXT J 
5770 ! 

5780 Totcor-Totcor-Conlen ! CONVERTS TOTAL CORRECTION FACTOR FROM <in> TO < 
mm) 

5790 * 

5800 ! 

5310 RETURN ! RETURNS TO Data acq OR TaKe_data ROUTINE 
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5820 f 

5830 

5840 ! 

5850 Crack.length: ! CALCULATES CRACK LENGTH USING HUDAK (OTMSRDC PROG) [PROF. J 
. JOYCE'S EQUATION] 

5860 ! NOTE: THIS ROUTINE USES COD COMPLIANCE METHOD VICE DEFLECTION COMPLIA 

NCE . 

5370 FOR I- I TO 6 
5880 READ Coef<I> 

5890 NEXT I 

5900 Comol-1/Stiff 

5910 C-2*H-Bef f -t-Compl /Hs/ 1 000 

5920 C-1 /<SGR(C>+1 > 

5930 G<1)-C 

5940 FOR 1-2 TO 5 

5950 G< I>-G< I- 1 >-G< T > 

5960 NEXT I 
5970 F*Coef < 1 ) 

5980 FOR 1-1 TO 5 

5990 F-F+Coef <1 + 1 )»G( I) 

6000 NEXT I 
6010 Crclcl-F-v4 
6020 RESTORE 

6030 DATA 0.998265.-3.81662,-1.80596,32.31033.-44.15665.-52.67876 
6040 ! 

6050 RETURN ! RETURNS TO Initial_a OR Main_prog 

6060 ! 

6070 ! 

6080 ! 



6090 

6100 

6T10 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

TAGE 



Conv_data: ! INPUTS CONVERSION DATA FOR MEASURING DEVICES 

PRINT CHRS( 12) ! CLEARS CRT 

PRINT "INPUT: LOAD RANGE. STRAIN RANGE. STROKE RANGE" 

PRINT " (ENTER IN PERCENT ie. 10.20,50,100)" 

INPUT Lr.Cr.Sr 
PRINT "" 

PRINT "LOAD RANGE - ”;Lr 

PRINT “STRAIN (COD) RANGE - ":Cr 

PRINT "STROKE (DEFLECTION) RANGE - ";Sr 

Lmajt-20000 ! MAX LOAD RANGE IN (lbs). 

Cmax- .200 ! MAX COD RANGE IN (in). 

Smax-5.0 ! MAX STROKE RANGE IN (in). 

V-10 ! MAX VOLTAGE IN (volts). 

* 

! MAXIMUM OPERATING POINTS ARE: 

01max-Lmax*Lr/1 00 
Ocmax -Cmax »C r / 1 0 0 
Qsmax-Smax-Sr/1 00 
* 



! CONVERSION CONSTANTS ARE 
Sle( 1 )-Olmax/V 
Sle(2)-0cmax/V 
Sie(3)-0smax/V 
! CONVERTING TO METRIC: 
Conlen-25 . 4 
Conib-4. 448221615 



(ENGLISH UNITS) 



! 25.4mm - 1 in 
! 4.448221615N - 



1 lb 



(HP PROGRAMING BOOK) 



! Sled. 2. 3) IS USED TO DENOTE CONVERSION OF VOLTAGE TO ENGLISH UNITS. 
SKI) -Sled )»Conlb*(-1 .0) ! COMPRESSIVE LOAD IS NEGATIVE VOTLAGE 

SI (2)-Sie(2)-Conlen 

SK3)-Sle(3)*ConlenM-1 .0) ! COMPRESSIVE STROKE SUPPLIES A NEGATIVE VOL 
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6 4 00 
0*4 1 U 

5420 

8430 
• • 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

66)0 

6620 

6530 

6640 

6650 

6660 

6670 

6630 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6320 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6980 

6990 






PRINT 

PRINT 

PRINT 



"YOUR ENTRIES HAVE PRODUCED THE FOLLOWING 



CONVERSION DATA 



PRINT "FOR LOAD: " :Si ( I ): "NEWTON/ VOLT (";$le< t > :"Lbs/V>" 
PRINT "" 

PRINT "FOR COD: " ;S1 (2) ;"mm/V <" :SLe<2> : "m/V>” 

PRINT 

PRINT "FOR DEFLECTION: M ;$l<3> ;"ran»/V <";Sie<3> :"in/V>" 
PRINT *•*• 

PRINT "WOULD YOU LIKE A HARD COPY?" 

INPUT HS 

IF H$-"N" OR HS-"NO" THEN 6630 

t 

PRINTER IS 706 ! TOGGLES PRINTER 

PRINT "" 

PRINT " CONVERSION DATA" 

PRINT "FOR L0AD:”;S1 < 1 ) ;"N/V <" ;SIe( 1 > :"Lbs/V)” 

PRINT "FOR C0D:":S1<2) ;"ram/V <":SIe<2> :"in/V>” 

PRINT "FOR DEFLECTI0N:”;SI<3) :"mm/V (" :Sle<3) ;"in/V>" 

PRINTER IS 1 ! TOGGLES PRINTER 

* 

PRINT CHRS(12) 

t 

RETURN ! RETURNS TO Ic-KEY SELECTION 



Slope: ! LINEAR LEAST SQUARES FIT 

PRINT CHRS< 12) ! CLEAR CRT 

! INITIALIZE VARIABLES 
P2-0 ! USED AS A COUNTER 

Suinx-0 ! SUM OF X VALUES 

Sumx2*0 ! SUM OF X SQUARED VALUES 

Siuny “0 * SUM OF Y VALUES 

Sumy2-0 ! SUM OF Y SQUARED VALUES 

Sujnxy-0 * SUM OF THE PRODUCT 

FOR I-K1 TO K3 ! Kt-INDEX OF FIRST. K3-INDEX OF LAST DATA POINT USED. 
! TYPE OF PLOT 

IF C2-1 THEN 6840 ! LOAD vs COD 

IF C2-2 THEN 6880 ! LOAD vs DEFL 

IF C2-3 THEN 6920 ! J vs DELTA_a 

I 

X“Cod< I) 

Y”Lod( I) 

GOTO 6990 

» 

X-DefKI) 

Y-Lod<I) 

GOTO 6990 

» 

X-Delta_a<I) 

Y»J9< I) 

» 

! CHECK IF DATA POINT IS WITHIN BOUNDS 

IF Del ta_a( IXJ9< I)/(2*Fs)* . 1 5 THEN 7050 * CHECKS LOWER LIMIT 

IF Delta_a(I)>J9<I)/<2-Fs)>1 .5 THEN 7050 ! CHECKS UPPER LIMIT 

* 

Sumx "Sumx+X 



77 



7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7130 

7200 

7210 

7220 

7230 



Sumy -Sumy ♦ f 
Sumx2-Sumx2 + X-X 



Sumy2-Sumy2 + Y*Y 
Sumxy-Sumxy+X-Y 

GOTO 7060 ! GET NEXT DATA POINT 

P2-P2+1 ! INCREMENTS COUNTER EACH TIME A POINT IS ELIMINATED. 



NEXT I 

K3-K3-P2 ? USED TO CALCULATE NUMBER OF DATA POINTS UITHIN BOUNDS. 

K4-K3-K1+1 ! K4-T0TAL * DATA POINTS USED IN CALC 

Oal -Sumx2-Sumx-*Sumx/K4 



Qa2-Sumxy-Sumx-Sumy /K4 
Qa3 -Sumy 2 -Sumy -Sumy /K 4 
; 



M-Qa2/Qa1 ! SLOPE 

» 

IF M<0 THEN 7430 * IF SLOPE < 0, BAD DATA 

! FROM DTNSRDC PROGRAM 
QI*<K4-Sumx2-Sumx'2)/(K4«(K4-| >) 
Q2-<K4«$umy2-Sumy'2>/<K4»<K4-l ) ) 

Q3*<K4-1 )»<02-M»M»01 >/<K4-2) 

Q4-SQR <Q3/<Q1 -<K4-1 ) ) ) 

» 

B-<Sumy-M-Sumx)/K4 ! INTERCEPT 



7240 Correl-M-SQR<Ql/Q2> 

7250 ! 

7260 Qa4-SGR<ABS< <Sumy2-B»Sumy-M»Suraxy >/<K4-2> ) ) 

7270 Sigmab-Qa4/SGR<ABS<Qa1 ) ) 

7280 Sigmam-Si.gn»aD»SQR<Sumx2/K4) 

7290 Stxff-M ! STIFF - SLOPE 

7300 Compl- 1 /M ! COD COMPLIANCE - 1 /SLOPE 

7310 ! 

7320 ! STUDENT t TEST 

7330 U4-Q4 

7340 T9- 1.96 

7350 IF K4> 1 20 THEN 7410 ! RETURNS 

7360 T9-2.1 

7370 IF K4<20 THEN 7410 ! RETURNS 

7380 N9-K4-2 

7390 T9-N9/(-.617109+N9-.5l0239) * 

7400 ! 

7410 RETURN 
7420 ! 

7430 PRINT "POOR DATA— SLOPE < 0" 

7440 M-333333 . 3 

7450 Correl-0 
7460 ! 

7470 RETURN ! RETURNS TO Inltial_a, Main_prog or R_curve ROUTINE (WHICHE 

VER CALLED IT) 

7480 ! 

7490 ! 

7500 * 

7510 Data.File: ! ESTABLISHES DATA FILE 

7520 ! THE FILE YOU CREATE CAN BE READ BY THE "J_DATA” PROGRAM 

7530 ! 

7540 ! NOTE: ANY CHANGES MADE TO THIS ROUTINE MAY REQUIRE SIMILAR CHANGES TO " 

J_DATA" PROGRAM DUE 

7550 ! TO DISK FORMATTING CONSULT HP MANUAL FOR ASSISTANCE. 

7560 ! 

7570 PRINT CHRS< 12) ! CLAERS CRT 
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7580 


PRINT "THIS ROUTINE HILL CREATE A DATA FILE GN DISK FOR 


ThE TEST." 


7590 


PRINT •••• 




7600 


PRINT "INSERT AN INITIALIZED OATA DISK AND ENTER FILE 


NAME (SAMPLE 


D IS 


SUGGESTED) . " 




7610 


INPUT NameS 




7620 


PRINT "" 




7630 


PRINT "IS THIS A NEW FILE?" 




7640 


INPUT OS 




7650 


IF Q$»"N" OR Q$-”NO" THEN 7680 




7660 


J 




7670 


CREATE BDAT NajneS.6.40000 




7630 


ASSIGN ^Fi.ie TO NameS 




7690 


j 




7700 


PRINT "" 




7710 


PRINT "FOR RECORD PURPOSES THE FOLLOWING DATA IS ALSO 


RECORDED:" 


? 720 


PRINT ”" 




7730 


PRINT "ENTER THE DATE (MN-DAY- YEAR) 




7740 


INPUT DaraS(l) . 




7750 


PRINT " M 




7760 


PRINT "ENTER THE TIME <HH:MM> 




7770 


INPUT DataS(2) 




7780 


PRINT 




7730 


PRINT ’•ENTER THE OPERATOR'S NAME :** 




7800 


PRINT " (UP TO 10 CHARACTERS)" 




7810 


INPUT DataS<3) 




7820 


PRINT "" 




7330 


PRINT "ENTER SPECIMEN ID:" 




7840 


PRINT " (UP TO 10 CHARACTERS)” 




7850 


INPUT DataS(4) 




7860 


! WRITE ON DISK 




7870 


FOR 1-1 TO 4 




7880 


OUTPUT «File;DataS( I)t 1 jlOl 




7890 


NEXT I 




7900 


• 




7910 


BEEP 




7920 


PRINT 




7930 


PRINT!- "DATA FILE HAS BEEN CREATED." 




7940 


; 




7950 


j 




'960 


RETURN ! RETURNS TO ic-KEY SELECTION 




7970 


* 




7980 






7990 


t 




8000 


Stor data: ! WRITES DATA ON DISK 




8010 


» . 




8020 


1 THE ITEMS STORED BY THIS ROUTINE CAN BE 




8030 


! RECALLED USING THE ”J DATA" PROGRAM 




8040 


! CARE MUST BE TAKEN IF ANY CHANGES ARE MADE TO THIS ROUTINE SUCH THAT : 


MILAR CHAGES ARE 




8050 


! MADE TO THE "J DATA" PROGRAM. 




8060 


t 




3070 


IF K2>0 THEN 8160 * K2 IS UNLOADING COUNTER 




8080 


OUTPUT £Fi le ;A0<D2> ♦ AminO <D2) , AmaxO<D2) , Comp 10 < 02) ,CorrelO(D2> ,K 


8090 


» 




8100 


PRINT "" 




8110 


PRINT "" 




8120 


PRINT "INITIAL CRACK LENGTH (a0> INFORMATION FROM 


MEASUREMENT * 


02 






8130 


PRINT "IS NOW IN DISK FILE ";NameS 




8140 


i 
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:3 1 50 IF <2*0 'HEN 8220 ! SKIP IF THIS WAS ad MEASUREMENT 

3160 OUTPUT »F i le :K2 . A 1 < K2 ) .Del ta_a< K2 > , K . L . Co r rel . Comp l . Am m . Amax . Area9 < K2 

) . J9(K2) 

3170 OUTPUT SFile:Mdefl<<2) .Codiff .K0.K1 

8180 * 

3190 PRINT "" 

3200 PRINT “UNLOADING INFORMATION IS NOW IN DISK FILE ”;NameS 
8210 ! 

8220 BEEP 

3230 ! WAIT 3 ! ALLOWS A PAUSE SO MESSAGE IS PRESENTED ON THE CRT 

8240 ! 

3250 ! 

8260 RETURN ! RETURNS TO Initial.a OR Main_prog 
8270 ! 

3280 - 

3290 * 

3300 Take data: ! THIS ROUTINE IS USED TO TAKE DATA DURING THE ACTUAL TEST. 

3310 * 

3220 I-< ! INDEX 

3330 OUTPUT Dvm:”AI1 VT1" ! READS LOAD 



8340 ENTER Dvm:Lod(I) 

3350 Lod(I)-Lod(I)*SI<1) 
3360 ! 

8370 OUTPUT Dvm:”AI2 VTt“ 
8380 ENTER Dvm:Cod(I> 

8390 Cod(I)*Cod(I)-Sl<2) 
8400 ! 

8410 OUTPUT Dvm:”AI3 VTt“ 
8420 ENTER DvmiDefKI) 

8430 Defl<I)-Defl(I)-Sl(3) 
8440 * 

3450 GOSUB Cor rect_def 1 
NCE 



! READS COD 



! READS DEFLECTION 



! DETERMINE CORRECTION FACTOR FOR SYSTEM COMPLIA 



8460 ! 

3470 Defi(I)-Oefl<I>-Totcor 
8430 ! 

3490 ! PLOT THE DATA POINT 

8500 VIEWPORT 13 . 100-RATI0 , 1 0 .90 

3510 WINDOW X1(C2).X2<C2).Y1(C2),r2<C2> 

3520 IF C2-2 THEN 8550 

3530 PLOT Cod< I) ,Lod<I> 

3540 GOTO 8560 

8550 PLOT DefKI).Lod(I) 

8560 ! 

3570 ! 

8580 RETURN ! RETURNS TO Main_prog 

8590 ! 

3600 ! ••*•••*•*•***•»*••**•«»••***••*****•••••*•*•••• 

8610 ! 

8620 J_caic: ! CALCULATES J WITH CRACK GROWTH CORRECTION ERNST MODEL (USED BY D 

TNSRDC) 

3630 ! 

8640 E5-2 

8650 ! CALCULATE J FOR FIRST <Jst> UNLOADING 

8660 IF K2>1 THEN 8690 

8670 J6-J7+E5/(W-AKK2) )»Area9<K2>/Bmin 

8680 GOTO 8730 

3690 JS”J7+E5/< W-AKK2) )*(Area9(K2) - Area9(K2- 1 > >/Bmin 

3700 ! IF DELTA_a < 0.254mm DO NOT USE CRACK GROWTH CORRECTION 

8710 IF Del ta_a(K2) < .254 THEN 8730 
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8720 J6-J6* 

STM E813-81 > 



< l-1/(W-Ai(K2)>*<Al(K2)-Al(K2-1>>> » CORRECTS FOR CRACK GROWTH <A 



8730 
8740 
8750 
8760 
8770 
8780 
3790 
8800 
8810 
8820 
8830 
8840 
8850 
EDING: 
8S60 
8870 
8880 
8890 
8900 
8910 
8920 
8930 
8940 
8950 
8960 
8970 
8980 
8990 
9000 
OKE' 
9010 
9020 
9030 
9040 



J7-J6 

i 



=J6 



? STORES CURRENT J VALUE IN J9 ARRAY. 



! RETURNS TO Main _prog 



J9(K2>> 

♦ 

RETURN 

\ 

Main_prog: ! THIS IS THE MAIN PROGRAM FOR 



CLEARS CRT 



THE J-INTEGRAL TEST. 



PRINT CHR$( 12) 

PRINT ,,M 

PRINT "ENSURE YOU HAVE COMPLETED THE FOLLOWING ITEMS 



BEFORE PROCE 



1 . 
2 . 

3. 

4. 



ENTERED SPECIMEN DATA." 

ENTERED CONVERSION DATA." 

SET UP INTEACTIVE PLOTTER." 
OBTAINED INITIAL CRACK LENGTH.' 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT HS 
IF HS»”NO 
PRINT "" 

PRINT 
PRINT 
INPUT 

Lod5=* . 

* 

PRINT 
CONTROL. 

PRINT 

PRINT "ENSURE THIS IS DONE." 
PRINT "" 

PRINT "WHEN YOU ARE READY PRESS 



TS MACHINE. 



ARE YOU READY?” 

OR HS«"N” THEN 10530 ! RETURNS TO WAIT FOR k-KEY SELECTION 



"WHAT IS YOUR EXPECTED MAXIMUM LOAD IN NEWTONS?" 

Pmax 

1 *Pmax ! USED TO SIGNAL TURN AROUND FROM UNLOADING 

"THE TEST SHOULD START WITH NEAR ZERO <0> LOAD WITH THE MTS IN 'STR 



'ENTER' OR 'CONTINUE' KEY AND START THE M 



9050 

9060 

9070 

9080 

9090 

9100 

9110 



! CLEARS CRT 



INPUT HS 
PRINT CHR$< 12) 

; 

* 

! INITIALIZE VARIABLES 

Area=0 ! AREA UNDER LOAD vs DEFLECTION CURVE DTNSRDC <I8> 
Area9< 1 >=0 ! AREA BENEATH SPECIFIC LOD vs DEFL (DTNSRDC I9<K2>> 



9120 


K-K 


? 


CURRENT INDEX (AFTER 


OBTAINING 


INITIAL CRACK LENGTH) 


9130 


I=K 


t 


INDEX OF FIRST 


TEST 


DATA 


POINT 




9140 


KO-K 


t 


WHERE TEST 


DATA 


STARTS 






9150 


K1 =0 


f 


LOWER INDEX 


FOR 


SLOPE CALC 




9160 


K2-Q 


t 


A COUNTER 




*» 


<K2> 


NUMBER OF UNLOADINGS 


9170 


K3*0 


} 


UPPER INDEX 


FOR 


SLOPE CALC 




9180 


L= 1 


t 


A COUNTER 


DTNSRDC 


<L> 


«• 


" DATA POINTS TAKEN 


UNLOADING 
















9190 


Lodi =*0 


t 


TEMP LOAD 




• • 


(PI ) 


USED 


IN AREA CALC 


9200 


Lod2=0 


t 


TEMP LOAD 




• l 


<P2> 


• t 


M M »• 


9210 


Lod5=* . 1 


♦Pmax ! 




•• 


<P5> 


USED 


FOR TURN AROUND SIGNAL 


9220 


Lod^O 


t 


MAX LOAD 




• t 


<P9> 






9230 


Codl =*0 


; 


TEMP COD 




M 


(DO 






9240 


Cod2-Q 




TEMP COD 




«• 


(D2> 






9250 


Cod4=Q 


l 


TEMP COD 




• • 


(C4) 
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9260 

9270 

9280 

9290 

9300 

9310 

9320 

9330 

9340 

9350 

9360 

9370 

9380 

9390 

9^00 

9410 

9420 

9430 

9440 

9450 

9460 

9470 

9480 

9490 

9500 

ADING 

9510 

9520 

9530 

9540 

9550 

9560 

9570 

9580 

9590 

9600 

9610 

9620 

9630 

9640 

9650 

9660 

9670 

9680 

9690 

9700 

9710 

9720 

9730 

9740 

9750 

9760 

9770 

9780 

9790 

9800 

9810 

9820 

9830 



Cod9=0 
Codif f ~0 
Def 1 1 =0 
Def 12=0 
J6=0 
J7=0 
11=0 
FI =0 

i 

» 

PRINT ’••• 



! MAX COD " (D9) 

! COD DIFFERENCE " <D3> 

! TEMP DEFL " <D4) USED IN AREA CALC 

! TEMP DEFL " <D5) 

! JI FROM CURRENT UNLOAD 
! JI FROM PREVIOUS UNLOAD DTNSDRC <J7> 

! FLAG < 0=L0ADING . 1 -UNLOADING) 

! FLAG (Ft =999 ENDS TEST) 



IF FI =999 THEN GOSUB End_test 
; 



GOSUB Take_data 
Def I6=Def 1<I> 

Lod7=Lod< I) 

Lod2=Lod( I) 

Def 12=Def KI) 

Cod2=Cod( I) 

t 

! TEST CONDITIONS OF DATA POINT JUST TAKEN 

IF (K-K0X20 THEN 10430 ! CALC AREA a TAKE MORE DATA 

Codiff=Cod2-Codl 

IF Codif f <-3 . OE-4 THEN 9540 ! IF CODIFF IS <-) CHECK IF LOADING OR UNLO 



IF 11=0 THEN 
IF Cod2>Cod9 
GOTO 10340 
IF I1O0 THEN 
Lod9=Lod2 
Cod9=Cod2 
Def 19=Def 12 
11 = 1 

GOTO 10340 
IF L<5 THEN 1 
KI-K-L 
K3-K1+L 
PENUP 

j 

K2-K2+1 

i 

GOSUB Slope 
Stif f2=Stif f 



10430 ! 

THEN 9600 
; 

1 0340 I 
* 

; 

» 



0290 



IF LOADING CALC AREA & TAKE MORE DATA 
! IF TRUE CHECK IF GOOD UNLOADING 
DISREGARD UPPER 27. OF UNLOADING CURVE 
IF UNLOADING CHECK IF ON UPPER 27. OF CURVE 
IF LOADING MAKE CURRENT LOAD = MAXIMUM LOAD 
MAKE CURRENT COD = MAX CCD 
MAKE CURRENT DEFL = MAX DEFL 
BE UNLOADING . 

TEST IF UPPER 27. OF UNLOADING 

! BAD UNLOADING, SET 11=1 

! SETS LOWER INDEX FOR SLOPE CALCULATION. 

! SETS UPPER INDEX FOR SLOPE CALCULATION. 



! INCREMENT UNLOADING COUNTER 



MUST 



! IF GOOD UNLOADING 



GOSUB Crack_length ! CALC CRACK LENGTH 
A1 (K2)=Crcic 1 ! CURRENT CRACK LENGTH 

Delta_a(K2)=Al(K2)-A1 ! DETERMINE CRACK GROWTH 

! 

! NOTE: A1 = INITIAL CRACK LENGTH; ASSIGNED IN Imtial_a ROUTINE. 



Crclc 14=Crclc 1 

PRINT " UNLOADING *";K2;" " 

! PRINT 

PRINT "CRACK LENGTH ="; A1<K2> 

PRINT "DELTA_a =";Delta_a<K2> ; ' V 
PRINT "TOTAL a DATA POINTS = ";K 
PRINT "» DATA POINTS THIS UNLOADING =";L 
PRINT "CORRELATION =";Correl 
PRINT "COD COMPLIANCE =” ;Comp l ;"mm/N" 
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9840 * 

9850 Area9< K2) “Area 
9860 ! 

9870 * CALCULATE 95% CONFIDENCE INTERVAL 

9830 St iff “Stiff + T9*U4 

9890 GOSUB Crack_length 

9900 Amajc-Crckl 

9910 Stif f=Stif f-2»T9*U4 

9920 GOSUB Crack_length 

9930 Stiff=Stiff2 

9940 Amin=Crckl 

9950 PRINT 

9960 PRINT "95% CONFIDENCE INTERVAL ON CRACK LENGTH FROM ";Amin;"TO";Amax 
9970 ! 

9980 GOSUB J_calc ! CALC J (USES J-ERNST MODEL) 

9990 * 

10000 Mdef l(K2)“Defl(K) • MAX DEFLECTION 
10010 ! 

10020 ! PRINT " M 

10030 PRINT "AREA =";Area9(K2> ;"N»mm" 

10040 PRINT "J “”;J9(K2) :"KJ/m~2” 

10050 PRINT "LOAD LINE DEFLECTION -";Mdef I(K2) ;"mm" 

10060 PRINT "CHANGE IN COD >";Codiff ;"mm" 

10070 PRINT ”KO“";KO;" K1“";K1 

10080 ! 

10090 • 

10100 GOSUB Prt_resalts ? PRINTS RESULTS ON PRINTER 

10110 ! 

10120 ! DATA POINTS- CAN BE PRESENTED FOR REVIEW BY SIMPLY REMOVING THE REMARK C 

HARACTER <!). 

10130 PRINT "" 

10140 ! PRINT " PT * LOAD(N) DEFLECTION(mm) COD(mm)" 

10150 ! FOR I»(K1-1) TO K 

10160 ! PRINT USING 9170:1. Lod(I> ,Def 1(1) .Cod(I) 

10170 ! IMAGE 4D,4X,5D.D,6X,DD.6D.7X,DD.6D 

10180 ! NEXT I 

10190 ! WAIT 5 ! PAUSE TO ALLOW VIEWING 

10200 ! 

10210 ! 

10220 GOSUB Stor_data ? STORES INFORMATION ON DISK 

10230 ! 

10240 ! 

10250 IF Crck 14-A1 >0 THEN 10290 
10260 PRINT "" 

10270 PRINT " NO CRACK EXTENSION THIS POINT." 

10280 PRINT "" 

10290 11=0 

10300 L=1 

10310 ! IF L=1 THEN K1=K ! SETS LOWER INDEX FOR SLOPE' CALCULATION AT CURRENT 

K 

10320 GOTO 10360 

10330 IF Cod2<Cod1 THEN 10440 ! GET MORE DATA 

10340 IF Lod2>.98»Lod9 THEN 10440 ! DISREGARD UPPER 2% OF UNLOADING 4 GET MO 

RE DATA 

10350 L*L+1 ! INCREMENT DATA POINT COUNTER 

10360 IF Lod2>Lod9-Lod5 THEN 10440 ! IF TRUE GET MORE DATA 

10370 BEEP ! SENDS OUT BEEP 

10380 BEEP ! AS TURN AROUND 

10390 BEEP ! SIGNAL TO BEGIN RELOADING 

10400 GOTO 10440 ! RESUME TAKING AND TEST DATA 

10410 ! 

10420 ! CALCULATE AREA UNDER LOAD vs DEFLECTION CURVE 
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10430 

10440 

10450 

10460 

10470 

10480 

10490 

10500 

10510 

10520 

10530 

10540 

10550 

10560 

10570 

10580 

10590 

10600 

10610 

10620 

10630 

10640 

10650 

10660 

10670 

10680 

10690 

10700 

10710 

10720 

10730 

10740 

10750 

10760 

10770 

10780 

10790 

10800 

10810 

10820 

10830 

10840 

10850 

10860 

10870 

10880 

10890 

10900 

10910 

10920 

10930 

10940 

10950 

10960 

10970 

10980 

10990 

11000 

11010 

11020 



Area=Area+ (Def 12-Def 1 1 >*< Lod2+Lod1 >/2 
Lodi =Lod2 

K-K + 1 ! INCREMENT TOTAL * DATA POINTS COUNTER 

Codl *Cod2 
Def 1 1 »Def 12 

f 

GOTO 9380 ! LOOP BACK CHECK IF TEST IS 

! OUFR . IF NOT TAKF MORF DATA. 



RETURN ! NEVER REACH THIS BUT ROUTINE WOULD NOT BE VALID WITHOUT IT 

t 4*****#*#*«»««*«*«****»*«»*««»*«**«»«tt*##» 

Prt_results: ! PRINTS TEST INFORMATION WHILE PROGRAM IS RUNNING 

PRINTER IS 706 ! TOGGLES PRINTER 

PRINT "" 

PRINT" UNLOADING *”;K2;” . " 

PRINT MM 

PRINT "CRACK LENGTH =":AKK2) :"mm DELTA _a =";Delta_a<K2> ;"mm" 

PRINT “95Z CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM" ;Anun;" TO”;Amax 
PRINT "" 

PRINT "TOTAL * DATA POINTS =";K,* M * DATA POINTS THIS UNLOADING =";L 

PRINT "CORRELATION »”;Correi 

PRINT "COD COMPLIANCE =";Compl ;"mm/N" 

PRINT "" 

PRINT "AREA -":Area9<K2> ;"N*mm” 

PRINT "J =”;J9(K2) ;"KJ/nT2" 

PRINT "MAX LOAD LINE DEFLECTION -";Mdef KK2) ;"mn»" 

PRINT "CHANGE IN COD »";Codif f ;"mm” 

; 

PRINT "" 

PRINTER IS 1 ! TOGGLES PRINTER TO CRT 

\ 

RETURN ! RETURNS TO Main_prog 

t 

j 

Stor_arrays: ! STORES ARRAYS (Lod.Def 1 ,Cod> ON DISK 

t 

! THE ARRRAYS CAN BE RECALLED LATER USING 
! THE "J_DATA M PROGRAM. 

OUTPUT 3Fiie;Lod<*),END 
OUTPUT SFile;Def 1<*>,END 
OUTPUT 3File;Cod<*> ,END 

t 

BEEP 

f 

♦ 

RETURN ! RETURNS TO End_test 




End_test : ! ENDS TEST 



\ 

PRINT CHRS< 12) 

\ 

\ 

GOSUB Stor_arrays 
\ 
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11030 
1 1 040 
11050 
11060 
11070 
11080 
11090 
11100 
11110 
11120 
11130 
11140 
11150 
11160 
11170 
11180 
11190 
11200 
11210 
1 1220 
11230 
11240 
11250 
IT 

11260 
1 1270 
11280 
11290 
11300 
11310 
11320 
11330 
11340 
11350 
11360 
11370 
11380 
11390 
11400 
11410 
11420 
11430 
11440 
11450 
11460 
11470 
11480 
11490 
11500 
11510 
11520 
11530 
11540 
11550 
I 1560 
11570 
11580 
11590 
11600 
11610 



PRINT "DATA HAS BEEN STORED ON DISK IN FILE ";Naroe5 

t 

GOSUB R_cu.rve 
» 

PRINT 

PRINT " TEST IS OVER!!” 

PRINT "" 

PRINT 
PRINT "" 

PRINT "" 

PRINT MM 

PRINT ” TEST IS OVER!” 

PRINT 
PRINT 
PRINT "" 

PRINT "" 

PRINT 

PRINT ” TEST IS OVER!” 



GOTO 12590 ! GO TO END STATEMENT 

» 

RETURN ! NEVER REACH THIS POINT BUT ROUTINE WOULD NOT BE VALID WITHOUT 



f -m- ■*■■*■**■■*■■* -m- ■*■*■*■*■ -tut ■*■*■•*■■*■■*■■**■ -m- 

Recurve: ! J vs DELTA_a DATA IS PRESENTED FOR PLOTTING 

f 

PRINT "” 

j 

! ALL <*) STATEMENTS WILL ALLOW DIRECT INPUT OF Delta_a(I) vs J9<I> DATA 
! YOU MUST REMOVE [!<»>] FROM BEGINNING OF EACH LINE TO ACTIVATE 
!(-) PRINT “INPUT THE * OF UNLOADINGS” 

!(*) INPUT K2 

! (») PRINT "ENTER DELTA_A< I> , J9< I)" 

PRINT "UNLOADING # DELTA_a(mm) J(kJ/nT2)" 

FOR 1-1 TO K2 

!(») INPUT Del ta_a( I) , J9< I) 

PRINT USING 1 1 420 : I , Del ta_a< I) , J9< I) 

IMAGE 6X,3D,4X,DD.6D,5X,6D.D 
NEXT I 
PRINT "” 

PRINT "DO YOU WANT A HARD COPY?” 

INPUT H£ 

IF HS-"N" OR H$-”NQ" THEN 11680 

; 

PRINTER IS 706 ! TOGGLES TO THERMAL PRINTER 

PRINT "UNLOADING * DELTA_a(mm) J < k J/nT2>" 

FOR 1-1 TO K2 

PRINT USING 1 1420;I,Delta_a<I) ,J9(I) 

NEXT I 

PRINTER IS 1 ! RETURNS TO CRT AS PRINTER 

» 

!<-> PRINT "ARE ALL OKAY?" 

! <*) INPUT GS 

!(») IF GS-”Y" THEN 11582 

!<*) PRINT "INPUT I (INDEX) OF BAD *" 

! (*> INPUT LI 

!(*) PRINT "INPUT VALUES OF DELTA_A( I) AND J9(I>" 
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11620 
1 ! 630 
11640 
11650 
11660 
A 

11670 
11680 
11690 
11700 
1 1710 
11720 
11730 
' OR ' 
11740 
11750 
11760 
1 1 770 
11780 
11790 
11800 
11810 
11820 
11830 
11840 
11850 
11860 
11870 
11880 
11890 
11900 
11910 
11920 
11930 
1 1940 
11950 
11960 
11970 
11980 
11990 
12000 
12010 
12020 
12030 
12040 
12050 
12060 
12070 
12080 
12090 
12100 
12110 
12120 
12130 
12140 
12150 
12160 
12170 
12180 
12190 



! ( *) INPUT De I ta_a< L 1 > , J9 < L i > 

!<-> PRINT "ARE YOU OK" 

! (»> INPUT EeeS 

!<»> IF Eee$-"N” THEN GOTO 11574 

1 END OF <»> CHANGES REQUIRED FOR DIRECT INPUT OF DELTA_A< I) vs J9(I> DAT 



PRINT "" 

PRINT "DO YOU WANT A PLOT OF J vs DELTA_a?" 

INPUT RS 

IF RS="N” OR R$-"NO" THEN 12550 
PRINT "" 

PRINT "WHEN THE PLOTTER IS READY (NEW GRAPH PAPER INSTALLED) PRESS 
CONTINUE' KEY" 

INPUT CS 



GOSUB Setup _p!otter 
; 

! PLOT THE DATA THE DATA POINTS 
VIEWPORT 13,100-RATIO, 10,90 
WINDOW XI (C2> ,X2(C2> , Y1 <C2> , Y2(C2> 
FOR 1-1 TO K2 



X-Del ta_a< I) 

Y-J9(I> 

MOVE X-X3-.005, Y+Y3-.005 1 DRAWS A 

DRAW X+X3- . 005 , Y- Y3* . 005 ! SMALL X 

MOVE X+X3-.005, Y+Y3-.005 * FOR DATA 

DRAW X-X3- .005, Y-Y3-.005 1 POINT 

NEXT I 
{ 

! PLOT THE BLUNTING LINE 
VIEWPORT 13,100-RATIO, 10,90 
WINDOW XI < C2 > , X2 < C2 > , Y 1 < C2 > . Y2 < C2 > 

Bly=Y2(C2>-Y3/8 ! TAKES LINE TO 7/8 OF Ymax 

Blx-Bly/(2»Fs>-1000 ! DELTA_a - J/2*Fs 

MOVE 0,0 
DRAW Blx ,Bly 



Blminx-Blx+. 15 
MOVE .15,0 
LINE TYPE 4 
DRAW Blminx,Bly 
; 

Blmaxx-Blx+I .5 
MOVE 1.5,0 
LINE TYPE 4 
DRAW Blmaxx.BIy 
PENUP 



! DRAWS LOWER BOUNDS 



1 DRAWS UPPER BOUNDS 



» 

K1-1 * LOWER INDEX FOR LEAST SQUARES FIT 

K3-K2 ! UPPER INDEX 

GOSUB Slope 



! PLOT THE LINE OBTAINED 
VIEWPORT 13.100-RATI0, 10,90 
WINDOW X1(C2) , X2(C2> ,Y1 <C2> . Y2(C2) 

Xplt=X2<C2>-X3/8 ! TAKES LINE TO 7/8 OF Xmax 

Yplt-Xplt-M+B 

MOVE 0 , B 

LINE TYPE 6 



'ENTER 



86 



12200 DRAW Xplt.Yplt 

12210 PENUP 

12220 PRINT "" 

12230 PRINT "ONLY DATA WITHIN BOUNDS IS USED" 

12240 PRINT "FOR LEAST SQUARES FIT." 

12250 PRINT •••• 

12260 PRINT "SLOPE’*" ;M:" INTERCEPT-" ;B 
12270 PRINT "CORREL-" :Correl 
12280 PRINT "NUMBER OF DATA POINTS USED:";K4 
12290 ! 

12300 PRINTER IS 706 ! TOGGLES PRINTER 

12310 PRINT'"" 

12320 PRINT "" 

12330 PRINT "DATA FROM J vs DELTA_a PLOT” 

12340 PRINT "" 

12350 PRINT "SLOPE -":M;" Y-INTERCEPT -";B 

12360 PRINT "CORRELATION -”;Correl;" FR0M";K4;" DATA POINTS." 

12370 PRINTER IS 1 ! TOGGLES PRINTER 

12380 ! 

12390 ! FIND THE INTERSECTION. HENCE J(1c>. 

12400 ! SOLVE BOTH EQUATIONS SIMULTANEOUSLY 

12410 ! J - M*(DELTA_a)+B - (DELTA_a)»(2*Fs> HENCE: 

12420 J1c-(B/M)/< 1 /M-< 1 000/ <2»Fs> ) ) 

12430 PRINT "" 

12440 PRINT " Jlc =";J1c 

12450 PRINTER IS 706 ! TOGGLES PRINTER 

12460 PRINT "" 

12470 PRINT " Jlc -";J1c 

12480 PRINT "" 

12490 PRINT "" 

12500 PRINT " TEST IS OVER!” 

12510 PRINTER IS 1 ! TOGGLES PRINTER 

12520 ! 

12530 ! GOTO 12100 ! GO TO END STATEMENT 

12540 ! 

12550 RETURN ! RETURNS TO End_test ROUTINE 

12560 ! 

1 257 0 * 

12580 I 

12590 END ! END OF TEST AND PROGRAM 

12600 ! 

1 26 10 ! 

12620 ! 

12630 SUB Stop_test ! STOPS TEST 

12640 ! Ic9-KEY ACTS AS PROGRAM INTERRUPT 

12650 ! Ic9-KEY IS NOT LABELED. BUT IS ASSIGNED A HIGHER PRIORITY TI INTERRUPT P 

ROGRAM 
12660 ! 

12670 COM /J_intgl/ FI 

12680 FI =999 ! SETS STOP FLAG 

12690 ! 

12700 SUBEND 
1 27 10 ! 
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APPENDIX F 



J DATA LISTING 



i o 



j_data program 



20 

30 

AO 

50 

60 

70 

80 

90 

too 
tto 
120 
130 
1 AO 
ISO 
160 
170 
180 
190 
200 
210 
220 
230 
2A0 
250 
260 
270 
280 
290 
300 
310 
320 
330 
3A0 
350 
360 
370 
380 
390 ' 
A00 
410 
420 
430 
440 
450 
460 
470 
480 
490 

J9(K2> 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 



THIS PROGRAM WILL ALLOW YOU TO READ A DATA 
FILE CREATED BY ”J_INTGL" . YOU MUST KNOW 
HOW MANY UNLOADINGS WERE RECORDED DURING 
THE TEST TO USE THIS PROGRAM. 

PROGRAM WAS WRITTEN BY: 

LCDR W.K. TRITCHLER 
VERSION 1.0 DATE 11/83 



DIM DataSd : 4 » C 101 

DIM Lod<8000> .Def K8000) ,Cod(8000) 

DIM AI< 100) .Area9d00> .Del ta_a< 1 00) 

DIM Mdef 1 < 1 00 ) , J9< 1 00) 

t 

PRINT CHRS<12) ! CLEARS CRT 

j 

DUMP DEVICE IS 706 

PRINT -WHAT DATA FILE DO YOU WISH TO READ?" 

INPUT NameS 

ASSIGN ?F lie TO NameS 

FOR 1-1 TO A 

ENTER ?File:Data$(I)C 1:101 
NEXT I 
PRINT "" 

PRINT "J- INTEGRAL TEST OF ";Data$<4> ;"BY: ";DataS<3> 

PRINT ”" 

PRINT "DATE: ~;DataS< 1 ) ;" TIME: “;DataS<2) 

ENTER §Fi le;Crclf 1 , Amin. Amax .Comp L .Correl ,K 
PRINT ”” 

PRINT — . 

PRINT "~ 

PRINT "aO-” ;Crck 1 ."WITH 95 X CONDIFENCE INTERVAL FROM";Anun;” TO":Amax 
PRINT "" 

PRINT "FROM COMPLIANCE - ";Compl 
PRINT 

PRINT "CORRELATION-” ; Correl ; "FROM” ;K; "DATA POINTS" 

PRINT "" 

PRINT "DO YOU WANT A HARD COPY OF THIS INFORMATION?" 

INPUT GS 

IF G$»"N” OR G$-"NO" THEN A50 

DUMP ALPHA 

PRINT 

PRINT "INPUT THE NUMBER OF UNLOADINGS THAT WERE RECORDED DURING THE TEST." 
INPUT Cl ! WILL BE USED AS COUNTER 

FOR I- I TO Ct 

ENTER $File;K2,AKK2> ,Delta_a<K2) , K,L, Correl .Compl .Anun.Amax ,Area9(K2) , 

ENTER SFile:Mdefl<K2>,Codlff.K0.K1 
GOSUB Prt_resul ts ! DISPLAYS INFO 
PRINT "" 

PRINT "DO YOU WANT A HARD COPY?" 

INPUT 0$ 

IF QS-"N" OR QS-”NO" THEN 580 
PRINTER IS 706 

GOSUB Prt results 1 WILL NOW PRINT INFO 
NEXT I 

PRINT "DO YOU WANT A LISTING OF THE DATA POINTS?" 
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600 INPUT E3 

610 IF E3-"N" OR E3-"N0" THEN 630 

620 PRINT 

630 PRINT "PT * LOAD (N) DEFLECTION <mm> COD (mm)" 

640 ENTER *Fi le:Lod( »> .Def l ( *> ,Coa( -> 

650 FOR I- I TO K 

660 PRINT USING 670 : 1 .Lod< I) .Def 1 < I) ,Cod< I) 

670 IMAGE X,4D.3X,5D.D,8X,D.6D.3X.D.6D 

680 NEXT I 

690 PRINT "" 

700 PRINT "DO YOU WANT A HARD COPY?" 

710 INPUT HS 

720 IF HS-"N" OR HS-"N0" THEN 830 

730 PRINTER IS 706 

740 PRINT "DATA FROM aO F0R:";DataS(4) ;"8Y: ";Data5(3> 

750 PRINT "" 

760 PRINT "TESTED ON: ";DataS< 1 ) -."TIME: ";Data$(2> 

770 PRINT "" 

780 PRINT "PT * LOAD <N> DEFLECTION <mm) COD <mm>" 

790 FOR 1-1 TO K 

300 PRINT USING 670 : I .Lod< I) . Def i< I) .Cod< I) 

810 NEXT I 

820 PRINTER IS 1 

830 PRINT CHRS< 12) ! CLEARS CRT 

840 PRINT "DO YOU WISH TO PLOT ANY DATA?" 

850 INPUT PS 

860 IF ES-"Y" OR ES-"YES" THEN 880 ! OTHERWISE MUST NOW READ DATA POINTS OFF 

DISKETTE 

870 ENTER ?Fi le:Lod< »> .Def I <• > ,Cod< ») 

880 IF PS»"Y" OR PS-"YES" THEN GOSUB Setup_pIotter 

890 PRINT "DO YOU WANT ANOTHER PLOT?" 

900 INPUT AS 

910 IF AS-"Y" OR AS- "YES" THEN GOSUB Setup_plotter 

920 PRINT "DO YOU WANT ANOTHER PLOT?" 

930 INPUT FS 

940 IF FS-"Y” OR FS-"YES" THEN GOSUB S*tup_plotter 
950 PRINT ” PROGRAM IS OVER!*!" 

960 PRINT "" 

970 PRINT "" 

980 PRINT "" 

990 PRINT "" 

1000 PRINT "" 

1010 PRINT " PROGRAM IS OVER!" 

1020 PRINT "" 

1030 PRINT "" 

1040 PRINT "" 

1050 PRINT "" 

1060 PRINT 

1070 PRINT " PROGRAM IS OVER!" 

1030 ! 

1090 GOTO 1290 
1100 ! 

1 1I o * 

1120 ! 

1130 Prt_results: ! DISPLAYS OR PRINTS UNLOADING INFORMATION 

1140 PRINT "" 

1150 PRINT "FROM UNLOADING *":K2:"REC0RDED";L:"DATA POINTS" 

1160 PRINT "COD COMPLIANCE -" :Comp l : "mm/N CORRELATION -"jCorrel 
1170 PRINT "CRACK LENGTH -":AI <K2) :"mm DELTA_A -";Dei ta_a(K2) :"mm" 

1180 PRINT "95X CONFIDENCE INTERVAL FROM": Amin: "mm TO" ; Amax : "mm" 
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; ! 90 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 



PRINT "AREA -";Area9<;<2> :"N-mm J 
PRINT "DATA POINT INDEXES FROM" :K1 ; "TO 
PRINT "" 

PRINTER IS 1 

» 

RETURN 



J9< K2> ; "KPa* 
” ;K : ” HERE U 



; 



Setup _p lot tor : ? SET UP PLOTTER 

! YOU CAN CHOOSE BETWEEN THREE DIFFERENT PLOTS 
PRINT CHRSU2) ! CLEARS CRT 
DIM TS[ 80 1 ,XS[80 1 , YSC80 1 
GCLEAR 
GINIT 

qhj 

PLOTTER IS 3.'*INTERNAL" ! PLOTS ON CRT 
PRINT "TYPE OF PLOT:" 

PRINT " 1. LOAD vs COD" 

PRINT " 2. LOAD vs DEFLECTION" 

PRINT ” 3. J vs CRACK EXTENSION" 

PRINT "ENTER THE * OF THE PLOT YOU DESIRE- 
INPUT C2 
PRINT "" 

PRINT "" 

IF C2-1 THEN 1470 
IF C2-2 THEN 1520 
IF C2-3 THEN t570 

TS«"PLOT of LOAD vs COD for "4DataS<4> 

PRINT TS 
XS-"COD <mm>" 

YS-"L0AD <N) " 

GOTO 1610 

TS--PLOT of LOAD vs DEFLECTION for "4DataS<4> 
PRINT TS 

XS-"DEFLECTION (mra)" 

Y$-"L0AD <N)" .. 

GOTO 1610 

TS-"PL0T OF J vs CRACK EXTENSION for "4DataS<4) 
PRINT TS 

X$-"DELTA_a (mm)” 

YS-" J <kJ/m*2>” 

PRINT "" 

PRINT "USING METRIC UNITS (mm.N.ic J/nT2.etc. ) 
PRINT "ENTER THE Xmm.Xma* .Ymin.Ymax" 

INPUT XI <C2) ,X2(C2) , Y1 <C2) , Y2(C2) 
X3-X2<C2)-X1(C2) 

Y3-Y2(C2)-YKC2) 

VIEWPORT 13.100-RATI0, 10, 90 
WINDOW XI <C2) ,X2<C2) , Y1 (C2) ,Y2(C2> 

PRINT 

PRINT "INPUT DELTA_X AND DEL TA_Y LABEL INTERVAL- 
PRINT "FOR THE TWO AXES:" 

INPUT X4,Y4 
PRINT "" 

PRINT "DO YOU WISH TO HAVE A GRID OVERLAY?" 

INPUT GS 

IF G$-"N" OR G$-"NO" THEN 1800 
PRINT "" 

PRINT "INPUT GRID LINES PER DELTA LABEL:" 
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O’) 3 



1 ’90 

1800 

1810 

1820 

1830 

1340 

1850 

18S0 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 



INPUT 80 

PRINT "PRESS 'ENTER' OR 'CONTINUE' KEY TO PLOT- 
INPUT HS ! DUMMY VARIABLE 

PRINT CHRS( 12) 1 CLEARS CRT 

IF G$-"N" OR G$-"NO" THEN 1970 
FOR I-X4/B0+XKC2) TO X2<C2> STEP X4/B0 
MOVE I,Y1(C2) 

DRAW I , Y2<C2> 

NEXT I 

FOR I“Y4/B0*Y! <C2> TO Y2<C2> STEP Y4/B0 
MOVE XI <C2) . I 
DRAW X2(C2),I 
NEXT I 
GOTO 1970 

IF GS«"Y" OR GS--YES" THEN 1970 
MOVE X2<C2) . Y1 <C2> 

DRAW X2<C2) ,Y2<C2) 

DRAW XI <C2> ,Y2<C2) 

CLIP ON 

AXES X4,Y4,X1<C2> ,YKC2) 

CLIP OFF 

VIEWPORT 13. 100-RATI0.5.90 
WINDOW XKC2),X2<C2) ,Y1 <C2> ,Y2<C2> 

LORG 4 
CsiZE 3 .6 

FOR X-Xt<C2> TO X2<C2> STEP X4 
CLIP OFF 
MOVE X.YKC2) 

LABEL USING -K";X ! PLOTS X-DIGITS 

NEXT X 

VIEWPORT 4.IOO-RATIO. 10,90 
WINDOW XI (C2> ,X2(C2) ,Y1 <C2> , Y2<C2) 

LORG 2 

FOR Y-YKC2) TO Y2(C2) STEP Y4 
CLIP OFF 
MOVE XI <C2) , Y 
LABEL USING "K";Y 
NEXT Y 

VIEWPORT 0,100-RATIO. 0. 100 
WINDOW XI <£2> ,X2<C2> , Yt <C2> ,Y2(C2) 

CSIZE 5 

LORG 6 ! PLOTS TITLE 

MOVE X2(C2)/2, Y2<C2) 

LABEL TS 

DEG 

LDIR 90 
CSIZE 4,. 7 

MOVE XI <C2> , Y2<C2>/2 
LABEL Y$ ! PLOTS Y-AXIS LABEL 

LORG 4 
LDIR 0 

MOVE X2(C2)/2,YKC2) 

LABEL XS ! PLOTS X-AXIS LABEL 

; 

PRINT "" 

PRINT — 

PRINT "" 

PRINT "" 

PRINT " IS THE PLOT (AXES AND LABELS) OK?” 

PRINT 
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2390 PRINT M NOTE: IF TOU ANSWER YES OR Y” 

2400 PRINT ” PLOT WILL BE SENT TO” 

2410 PRINT ” EXTERNAL PLOTTER” 

2420 INPUT FS 

2^30 IF FS-”N” OR FS-"NO" THEN GOSUB Setup .plotter 1 REDO PLOTTER SETUP 
2440 ! 

2450 GCLEAR l CLEARS PLOT FROM CRT 

2460 ! 

2^70 G0SU8 Ext_plotter 
2480 1 

2490 PENUP 
2500 MOVE 0,0 
25 1 0 f 

2520 GOSUB Plot.data 
2530 ! 

2540 RETURN 
2550 ! 

2560 ! 

2570 ! 

2580 Ext_plotter : ! DUMPS PLOTTER SETUP ON EXTERNAL HP-7225B PLOTTER 

2590 PLOTTER IS 705,”HPGL” 

2600 VIEWPORT 13. 100-RATIO, 10 .90 
2610 WINDOW X1(C2).X2(C2).Y1<C2).Y2<C2) 

2S20 IF G3-”N” OR GS-”N0“ THEN 2730 

2630 FOR I-X4/B0+XKC2) TO X2(C2> STEP X4/B0 

2640 MOVE I,Yt<C2) 

2650 DRAW I,T2(C2) 

2660 NEXT I 

2670 FOR I-Y4/80+Y1 (C2) TO Y2<C2) STEP Y4/B0 

2680 MOVE XI (C2) ,1 

2690 DRAW X2<C2),I 

2700 NEXT I 

2710 GOTO 2760 

2720 IF GS-”Y” OR G$-”YES“ THEN 2760 
2730 MOVE X2(C2),Y1<C2) 

2740 DRAW X2(C2) . Y2<C2) 

2750 DRAW XI <C2) . Y2<C2) 

2760 CLIP ON 

2770 AXES X4.Y4,XJ<C2).Y1<C2) 

2780 CLIP OFF 

2790 VIEWPORT 13, 100-RATI0, 5, 90 

2800 WINDOW X1(C2) ,X2<C2> ,Y1 <C2) ,Y2(C2) 

2810 LORG 4 
2820 CSIZE 3.. 6 

2830 FOR X-XKC2) TO X2(C2) STEP X4 
2840 CLIP OFF 

2850 MOVE X.YKC2) 

2860 LABEL USING ”K”;X 

2870 NEXT X 

2880 VIEWPORT 4, tOO-RATIO, 10 .90 

2890 WINDOW XI <C2> ,X2<C2> . Y1 <C2> , Y2(C2) 

2900 LORG 2 

2910 FOR Y-YKC2) TO Y2(C2> STEP Y4 
2920 CLIP OFF 

2930 MOVE XI (C2> , Y 
2940 LABEL USING ”K";Y 

2950 NEXT Y 

29&0 VIEWPORT 0. TOO-RATIO. 0. 100 

2970 WINDOW XI <C2) ,X2<C2) . Y1 (C2) , Y2(C2) 

2980 LORG 2 
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2990 CSIZE 5 
3000 LORG 6 

3010 MOVE X2<C2)/2.Y2<C2> 

3020 LABEL TS 

3030 DEG 

3040 LDIR 90 

3050 CSIZE 4 ,. 7 

3060 MOVE XI <C2) . Y2<C2)/2 

3070 LABEL YS 

3080 LORG 4 

3090 LDIR 0 

3100 MOVE X2(C2>/2, Yt <C2) 

3110 LABEL XS 

3120 ! 

3130 PENUP 
3140 MOVE 0.0 
3150 ! 

3160 RETURN 
3170 ! 

3180 ».*..*****.**.***« 

3190 ! 

3200 PIot_data: ! PLOTS TO DATA SELECTED 

3210 ! 

3220 VIEWPORT 13 . 100-RATI0. 1 0 .90 
3230 WINDOW XI (C2> ,X2<C2) . Y1 <C2> , Y2<C2> 

3240 ! 

3250 IF C2* I THEM 3290 
3260 IF C2-2 THEN 3340 
3270 IF C2-3 THEN 3390 
3280 ! 

3290 FOR I»1 TO K 

3300 PLOT CodC I) ,Lod<I) 

3310 NEXT I 
3320 GOTO 3430 
3330 ! 

3340 FOR I-f TO K 

3350 PLOT Def KI) ,Lod<I> 

3360 NEXT I 
3370 GOTO 3430 
3380 ! 

3390 FOR I- I TO Cl 

3400 PLOT Del ta_a<K2> . J<K2> 

3410 NEXT I 
3420 ! 

3430 PENUP 
3440 BEEP 

3450 RETURN f RETURNS TO SETUP_PLOTTER 

3460 ! 

3470 I——.—..——.—..——.————— 

3480 ! 

3490 END 
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APPENDIX G 



PRECRACK OUTPUT 



USING ENGLISH UNITS dn.lbs.psi> THE FOLLOWING INFORMATION 
HAS BEEN ENTERED FOR: HY80-5B 

WIDTH * .9469 NOTCH LENGTH - .444 

Bmax » .496 Bmin * .394 

HALF SPAN - 1.894 POISSON'S RATIO * .3 • 

YIELD STRENGTH - 86600 ULTIMATE TS - 100000 
ELASTIC MODULUS - 2.9E+7 FLOW STRESS • 93300 

CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) FOR: HY80-5B 
WIDTH - 24.05126 NOTCH LENGTH * 11.2776 

Bmax * 12.5984 Bmin * 10.0076 

YIELD STRENGTH - 597085.603738 ULTIMATE TS - 689475.293 
ELASTIC MODULUS » 1 .9994783497E+8 FLOW STRESS * 643280.448369 

BASED ON CURRENT INFORMATION FOR HY80-5B: 

CURRENT CRACK LENGTH « .444 in 

OPTIMUN CRACK LENGTH - .56314 in 

EXTENSION REQUIRED » .12414 in 

CRACK LENGTH AT 501' OF EXTENSION - .50607 in 

Pmax » 3272.42651647 Lbs Pmm - 163.621325823 Lbs 

Dele - 3108.80519064 Lbs Kmax - 32563.7466931 

DELTA_K - 78435.5593585 in*. 5 DELTA_KE » 40000 in*. 5 

NEW LOADS BASED ON DELTA.KE ARE: 

Pmax - 1668.84843723 Lbs Pmm - 83.4424218615 Lbs 
DELTA_P - 1585.40601537 Lbs Kmax - 42105.2631579 



CONVERSION DATA 
FOR LOAD : - 1 779.288646 N/V < 
FOR COD: .0508 mm/V ( .002 
FOR DEFLECTION: -12. 7 mm/V ( 
U - .07 

a/W - .732186117982 in CRACK 
LOAD < lbs) * COD < in> 



200 


.00256290169302 


300 


.00384435253953 


400 


.00512580338604 


500 


.00640725423255 


600 


.00768870507906 


700 


.80897015592557 


800 


.0102516067721 


900 


.0115330576186 


1000 


.0128145084651 


1100 


.0140959593116 


1200 


.0153774101581 


1300 


.0166588610046 


1400 


.0179403118511 


1500 


.0192217626976 


1600 


.0205032135441 


1700 


.0217846643907 


1800 


.0230661152372 


U - .075 




a/W - . 7t 3968763392 in CRACK 


LOAD (lbs) 


COD (in) 


200 


.00220863052208 


300 


. 0033129457831 2 


400 


.00441726104416 


500 


.00552157630521 



5 J 



400 Lbs/V) 

n/V) 

.5 in/V) 

LENGTH - .693307035117 in 
COD (Volts) 

1 .28145084651 

I . 92217626976 
2.562901693-32 
3.20362711627 
3.84435253953 
4.48507796278 
5.12580338604 
5.76652880929 
6.40725423255 
7.0479796558 
7.68870507906 
8.32943050231 
8.97015592557 
9.61088134832 
10.2516067721 
10.8923321953 

II. 5330576186 

LENGTH - .676057022056 in 
COD (Volts) 

1.10431526104 
1 .65647289156 
2.20863052208 
2.7607881526 

5 
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USING ENGLISH OMITS <in.Ibs.psii IKE FOLLOWING INFORMATION 
HAS 3EEN ENTERED FOR: HY80-53 

WIDTH - .3-569 NOTCH LENGTH ■ .444 

3max * .4S6 8mm » .394 

hALF SPAN * 1.394 POISSON'S RATIO * .3 
YIELD STRENGTH * 36600 ULTIMATE TS - 100000 
ELASTIC MODULUS - 2.9E+7 FLOW STRESS - 93300 

CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPai FOR: HY80-5B 
WIDTH - 24.05126 NOTCH LENGTH - 11.2776 

3 max * 12.5984 3mm ■ 10.0076 

YIELD STRENGTH - 597085.603733 ULTIMATE TS - 689475.293 
ELASTIC MODULUS - 1 .9994783497E+8 FLOW STRESS > 643280.448369 

BASED ON CURRENT INFORMATION FOR HY80-5B: 

CURRENT CRACK LENGTH * .444 m 

GPTIMUN CRACK LENGTH - .56314 m 

EXTENSION REQUIRED * 712414 m 

CRACK LENGTH AT SOX OF EXTENSION » .50607 m 

Pmax - 3272.42651647 Lbs Pmm » 163.621325323 Lbs 

Delp - 3103.80519064 Lbs Kmax - 82563.7466931 

DELTA_K - 78435.5593585 in'. 5 DELTA_KE - 40000 in*. 5 

NEW LOADS BASED ON DELTAJKE ARE: 

Pmax - 1668.84843723 Lbs Pmm * 83.4424218615 Lbs 
DELTA_P > 1585 .40601537 Lbs Kmax - 42105.2631579 



CONVERSION DATA 

FOR L0AD:-1 779.283646 N/V <-400 Lbs/V) 

FOR COD: .0508 mm /V < .002 in/V) 

FOR DEFLECTION :-1 2 . 7 mm /V <-.5 in/V) 

THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .477215195165 m FROM MEASUREMENT * 1 
WITH 95 '/.CONFIDENCE INTERVAL FROM .47792545433 TO .476507410981 in 
USING COD COMPLIANCE OF: 2 .841 0431 1 054E-6 m/Lb - 1 .629051 73472E-5 mm/N 
WITH A CORRELATION OF: .999771386472 FROM 106 DATA POINTS. 



THE FOLLOWING INFORMATION HAS 3EEN OBTAINED: 

CRACK LENGTH - .477198279382 in FROM MEASUREMENT * 2 
WITH 95XCONFIDENCE INTERVAL FROM .478064054454 TO .476336178981 in 
USING COD COMPLIANCE OF: 2 .840761 26094E-6 m/Lb ■ 1 .6303775021 3E-5 mm/N 
WITH A CORRELATION OF: .99959257929 FROM 124 DATA POINTS. 

BASED ON CURRENT INFORMATION FOR HY80-5B: 

CURRENT CRACK LENGTH ■ .477198279382 in 

OPTIMUN CRACK LENGTH - .56314 m 

EXTENSION REQUIRED - .0909417206182 in 

CRACK LENGTH AT 50% OF EXTENSION - .489470860309 in 

Pmax - 2854.63728597 Lbs Pmin • 142.731864298 Lbs 

Deip * 2711.90542167 Lbs Kmax - 80310.2753096 

DELTA_K > 76294.7615441 m*.5 DELTA_KE - 40000 in .5 

NEW LOADS BASED ON DELTA_KE ARE: 

Pmax * 1496.63606161 Lbs Pmm * 74.8313030804 Lbs 
DELTA P - 1421 .30425353 Lbs Kmax - *2105.2531579 
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*jL_jWu1G INFORMATION ~HS SEEN OBTAINED : 

CRACK.' LENGTH - .499315ini338 m FROM MEASUREMENT * 1 
WITH 357.CONFIDENCE INTERVAL FROM .500434044632 TO .498202259284 in 
USING COD COMPLIANCE OF: 3 .23884557474E-6 m/Lb • 1 .861 83031 428E-5 mm/M 
WITH A CORRELATION OF: .399477933275 FROM 34 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH * .499773214641 in FROM MEASUREMENT * 2 
WITH 957.C0NFIDENCE INTERVAL FROM .50177796077 TO .497787999797 in 
USING COD COMPLIANCE OF: 3 .24775972936E-6 m/Lb * 1 .8769892409E-5 mm/M 
WITH A CORRELATION OF: .999019973378 FROM 57 DATA POINTS. 

BASED ON CURRENT INFORMATION FOR HY80-58: 

CURRENT CRACK LENGTH * .499773214641 in 

OPTIMUN CRACK LENGTH - .56314 m 

EXTENSION REQUIRED - .0633667853593 in 

CRACK LENGTH AT 507. OF EXTENSION * .47818339268 in 

Pmax ■ 1034.73227996 Lbs Pmm ■ 51.736613998 Lbs 

Delp - 982.995665963 Lbs Kmax - 31462.7523142 

DEL T A _K - 29889.6146935 m".5 DELTA.KE - 29000 in*. 5 

NEW LOADS BASED ON DELTA_KE ARE: 

Pmax * 1003.9351936 Lbs Pmm • 50.1967596798 Lbs 
0ELTA_P - 953.738433917 Lbs Kraax * 30526.3157895 



THE FOLLOWING INFORMATION HAS BEEN 08TAINED: 

CRACK LENGTH - .523238325069 in FROM MEASUREMENT * 1 
WITH 952C0NFIDENCE INTERVAL FROM .52495282665 7 TO .521538t09085 in 
USING COD COMPLIANCE OF: 3 . 74604099481 E-6 in/Lb - 2 . 1 61 81 567648E-5 mm/N 
WITH A CORRELATION OF: .999183456903 FROM 61 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .519260644793 in FROM MEASUREMENT * 2 
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WITH 957.C0NFIDENCE INTERVAL FROM .521257634776 TO . 5 1 72330 1 U25 in 
USING COD COMPLIANCE OF: 3 .65545467572E-S in/Lb - 2 . 1 1 309694446E-5 mm/M 
WITH A CORRELATION OF: .998905179767 FROM 61 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH =■ .520164834896 in FROM MEASUREMENT * 3 
WITH 957.C0NFIDEMCE INTERVAL FROM .521 732160051 TO .51360945694 in 
USING COD COMPLIANCE OF: 3.67581279234E-6 in/Lb » 2 . 1 1 927866299E-5 mm/N 
WITH A CORRELATION OF: .999107091656 FROM 79 DATA POINTS. 

3ASED ON CURRENT INFORMATION FOR HY80-5B: 

CURRENT CRACK LENGTH - .520164834396 in 

OPTIMUN CRACK LENGTH ■ .56314 in 

EXTENSION REQUIRED - .0479751651041 in 

CRACK LENGTH AT 507. OF EXTENSION - .467987582552 in 

Pmax * 942.504620395 Lbs Pmin - 47.1252310197 Lbs 

Delp - 895.379389375 Lbs Kmax =■ 30831.6649108 

DELTA.K » 29290.0816653 in*. 5 DELTA.KE - 29000 in*. 5 

NEW LOADS BASED ON DELTA.KE ARE: 

Pmax * 933.170289649 Lbs Pmin - 46.6585144825 Lbs 
DELTA_P » 886.511775167 Lbs Kmax - 30526.3157895 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH =■ .557542180253 in FROM MEASUREMENT * 1 
WITH 957.CONFIDENCE INTERVAL FROM .559342565208 TO .555757596637 in 
USING COD COMPLIANCE OF: 4 .6522934609E-6 in/Lb - 2 .6377228564E-5 mm/N 
WITH A CORRELATION OF: .999026349545 FROM 60 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .556985314799 in FROM MEASUREMENT * 2 
WITH 95XC0NFIDENCE INTERVAL FROM .559603140017 TO .554400760935 in 
USING COD COMPLIANCE OF: 4 .63555555373E-6 in/Lb * 2.69227537518E-5 mm/N 
WITH A CORRELATION OF: .998286279816 FROM 51 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .55746232304 in FROM MEASUREMENT * 3 
WITH 95XCONFIDENCE INTERVAL FROM .558944733365 TO .555990641885 in 
USING COD COMPLIANCE OF: 4 .6498885651 2E-6 in/Lb - 2.68079074153E-5 mm/N 
WITH A CORRELATION OF: .999219985095 FROM 70 DATA POINTS. 

BASED ON CURRENT INFORMATION FOR HY80-5B: 

CURRENT CRACK LENGTH - .55746232304 in 

OPTIMUN CRACK LENGTH - .56814 in 

EXTENSION REQUIRED - .0106776769602 in 

CRACK LENGTH AT 507. OF EXTENSION - .44933883348 in 

Pmax - 784.95101444 Lbs Pmin - 39.247550722 Lbs 

Delp - 745.703463718 Lbs Kmax - 29594.1425953 

DELTA_K - 28114.4354655 in*. 5 DELTA.KE » 29000 in*. 5 

THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .556873246958 m FROM MEASUREMENT * 1 

WITH 957CONFIDENCE INTERVAL FROM .558630143876 TO .555131398279 in 

USING COD COMPLIANCE OF: 4.63219612139E-6 in/Lb - 2 . 6753248 1 877E-5 mm/N 
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WITH A CORRELATION OF: . 33S 1 741 71419 FROM 5*4 DATA POINTS. 

BASED ON CURRENT INFORMATION FOR HY30-58: 

CURRENT CRACK LENGTH * .556873246958 in 

CPTIMUN CRACK LENGTH * .56814 in 

EXTENSION REQUIRED - .0112667530417 m 

CRACK LENGTH AT 507. OF EXTENSION * .449633376521 m 

Pmax - 787.327495262 Lbs Pmin - 39.3663747631 Lbs 

D*?1p - 747.961 120499 Lbs Kmax * 29614.5202835 

DELTA.K « 28133. 30376S4 in*. 5 DELTA_KE *■ 23000 in*. 5 

THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH * .560000734225 m FROM MEASUREMENT * 1 
WITH 95 /.CONFIDENCE INTERVAL FROM .56106384425 TO .55893721925 in 
USING COD COMPLIANCE OF: 4 . 7270940231 8E-6 in/Lb - 2 . 71 807983847E-5 mm/N 
WITH A CORRELATION OF: .99933646638 FROM 101 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .562208107534 in FROM MEASUREMENT * 2 
WITH 957.CQNF IDENCE INTERVAL FROM .56343540602 TO .560988180658 in 
USING CCD COMPLIANCE OF: 4. 79552996551E-5 in/Lb ■ 2 . 76035025025E-5 mm/M 
WITH A CORRELATION OF: .99948124671 FROM 67 DATA POINTS. 

BASED ON CURRENT INFORMATION FOR HY80-58: 

CURRENT CRACK LENGTH - .562208107534 in 

OPTIMUN CRACK LENGTH - .56814 in 

EXTENSION REQUIRED * .00593139246594 in 

CRACK LENGTH AT 507. OF EXTENSION * .446965946233 m 

Pmax - 765.936366326 Lbs Pnun > 38.2968183162 Lbs 

Delp - 727.63954801 Lbs Kmax - 29428.899829 

DELTA_K * 27957.4548376 in'. 5 DELTA_KE - 29000 in'. 5 

THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH = .564151471472 in FROM MEASUREMENT * 1 
WITH 957.CONFIDENCE INTERVAL FROM .565674835605 TO .562639453926 in 
USING COD COMPLIANCE OF: 4 .85680462792E-6 in/Lb - 2 .801 12437379E-5 mm/M 
WITH A CORRELATION OF: .999133273453 FROM 72 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

CRACK LENGTH - .563361917225 in FROM MEASUREMENT * 2 
WITH 957.C0NFIDENCE INTERVAL FROM .569305393318 TO .567422811142 in 
USING COD COMPLIANCE OF: 4 .99294827909E-6 in/Lb - 2.86883916896E-5 mm/N 
WITH A CORRELATION OFi .99960 J697254 FROM 84 DATA POINTS. 

PRECRACKING IS COMPlEIED! 

FINAL CRACK LENGTH - .568361917225 
CORRELATION OF DATA * .399601697254 

REMOVE SPECIMEN FROM MTS. 
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APPENDIX H 



J INTGL OUTPUT 



THREE POINT BEND SPECIMEN 
TEST REPORT 

SPECIMEN I.D. W/gQ-sS'vS MATERIAL j // Q Q ORIENTATION T- S 

Data /£//>7/ft.3 Operator ^7^ < tc-M <- C^f Teat Teao ^2 / d 



TT 

Em ax 

_i_L 




t4 ff 

U - A * a ^4 b . 0. 503- 9 L - 4, 

\ 

Ha * /. 8 9 4* Qmax » o-49e, o' Bain * 4 i 

Ya * @ b hGC) y 2X “ 1 * Uta * JOO &Q & jQu**- Fa * 9^ 3 O <2 

Poiaaon'a ratio * ^ Elaatic Modulus * X » 9 </fl 

Consents: 

— H.y BP S/«ftC- tUsJZAJz EsLOM PyLMSiAXtJl hlt^x l 

ft«J /a //</f.? ' 3*- Tc-st /^//,i/e3 
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USING ENGLISH UNITS <jn.Ibs.psi> THE FOLLOWING INFORMATION 
HAS BEEN ENTERED FOR: HYSG-5B 

WIDTH * .9469 NOTCH LENGTH * .444 

Bmax " .496 Bmin * .394 

HALF SPAM » 1.394 POISSON'S RATIO * .3 
YIELD STRENGTH - 86600 ULTIMATE TS - 100000 
ELASTIC MODULUS - 2.9E+7 FLOW STRESS - 93300 

CONVERTING THIS DATA TO METRIC UNITS (mm.N.KPa) FOR: HY80-5B 

WIDTH - 24.05126 NOTCH LENGTH = 11 .2776 

Bmax » 12.5984 Bnun * 10.0076 

HALF SPAN - 48.1076 POISSON'S RATIO - .3 

YIELD STRENGTH - 5970S5 .603738 ULTIMATE TS - 639475.293 

ELASTIC MODULUS - 1 .9994783497E+8 FLOW STRESS = 643280.443369 



CONVERSION DATA 

FOR LOAD : -1 779.288646 N/V < 400 Lbs/V> 

FOR COD: .254 mm/V < .01 in/V> 

FOR DEFLECTION :-1 .27 mm/V < .05 m/V) 

GIVEN a- .568 inches THEM: 

Pmax- 895.290182804 Lbs FROM ASTM E813-81 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

INITIAL CRACK LENGTH - 14.3875332557 mm FROM MEASUREMENT * 1 

WITH 957. CONFIDENCE INTERVAL FROM 14.4079726404 TO 14.3671845862 mm 

USING COD COMPLIANCE OF: 2 .81 520064535E-5 mm/N 

UITH A CORRELATION OF: .999323843963 FROM 189 DATA POINTS. 



THE FOLLOWING INFORMATION HAS BEEN OBTAINED: 

INITIAL CRACK LENGTH - 14.4326755548 mm FROM MEASUREMENT a 2 

WITH 957. CONFIDENCE INTERVAL FROM 14.4442813846 TO 14.42109582 mm 

USING COD COMPLIANCE OF: 2.84829954475E-5 mm/N 

WITH A CORRELATION OF: .99984474937 FROM 134 DATA POINTS. 
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INITIAL CRACK LENGTH * 1 4 . 445571 2237 mm FROM MEASUREMENT * 3 

WITH 95% CONFIDENCE INTERVAL FROM 1 4.4565506791 TO 14.434615129 mm 

USING COD COMPLIANCE OF: 2 .8578451 4764E-5 mm/N 

WITH A CORRELATION OF: .999816522357 FROM 176 DATA POINTS. 



UNLOADING * 1 

CRACK LENGTH » 14.4487507815 mm DELTA _a - .016075226716 mm 

95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.4693514291 TO 14.4282322251 

TOTAL * DATA POINTS - 632 * DATA POINTS THIS UNLOADING - 128 

CORRELATION > .999529434569 

COD COMPLIANCE - 2 .8755529381 7E-5 mm/N 

AREA - 711.365042717 N*mm 

J - 14.8049806501 KJ/m‘2 

MAX LOAD LINE DEFLECTION * .24034834255 mm 

CHANGE IN COD - .00109982 mm 



UNLOADING a 2 

CRACK LENGTH - 14.4711513219 mm DELTA _a - .0384757670308 mm 

95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.4927078135 TO 14.4496847244 

TOTAL * DATA POINTS * 947 # DATA POINTS THIS UNLOADING - 132 

CORRELATION - .9994671 1 1 958 

COD COMPLIANCE - 2. 8930791 3363E-5 mm/N 

AREA - 1693.69098454 N*mm 
J « 35.2970225265 KJ/m‘2 

MAX LOAD LINE DEFLECTION » .382347556072 mm 
CHANGE IN COD - .0008382 mm 



UNLOADING * 3 

CRACK LENGTH - 14.4875164202 mm DELTA _a - .0548408654168 mm 

95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.5143346851 TO 14.460837178 

TOTAL * DATA POINTS -1300 * DATA POINTS THIS UNLOADING - 157 

CORRELATION - .999017363069 

COD COMPLIANCE » 2.90942709989E-5 mm/N 

AREA = 3033.0003923 N*mm 
J - 63.28381022 KJ/m'2 

MAX LOAD LINE DEFLECTION - .543433314464 mm 
CHANGE IN COD - .0009398 mm 



UNLOADING #4 

CRACK LENGTH * 14.5414951598 mm DELTA _a - .108819604992 mm 

95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.5601316376 TO 14.5223259715 

TOTAL * DATA POINTS = 1670 f DATA POINTS THIS UNLOADING - 109 

CORRELATION * .999675611283 

COD COMPLIANCE - 2.944429881 18E-5 mm/N 

AREA - 5208.32927275 N-mm 
J * 108.998439815 KJ/m'2 

MAX LOAD LINE DEFLECTION * .782918890293 mm 
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CHANGE IN COD - .0007874 mm 



UNLOADING » 5 

CRACK LENGTH - 1 4 . 7407494376 mm DELTA _a - .308073882773 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.7587196422 TO 14. 722841991 1 

TOTAL * DATA POINTS » 2060 * DATA POINTS THIS UNLOADING - 95 

CORRELATION - .999731105207 

COD COMPLIANCE - 3 . 1 0248861 726E-5 mm/N 

AREA - 82t6. 7400787 N»mm 
J = 169.858689852 KJ/nT2 

MAX LOAD LINE DEFLECTION - 1.10455056629 mm 
CHANGE IN COD - .0008636 mm 



UNLOADING * 6 

CRACK LENGTH - 14.8152586049 mm DELTA _a - .382583050059 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.8481840789 TO. ,14. 782543466 

TOTAL * DATA POINTS - 2280 * DATA POINTS THIS UNLOADING - 138 

CORRELATION = .998642018089 

COD COMPLIANCE - 3 . 1 771 921 9461 E-5 mm/N 

AREA - 9001.97373669 N-mm 
J - 185.342172083 KJ/nT2 

MAX LOAD LINE DEFLECTION - 1.18989063161 mm 
CHANGE IN COD - .0009906 mm 



UNLOADING * 1 

CRACK LENGTH » 14.8695206824 mm DELTA _a - .436845127588 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 14.3849566722 TO 14.8541311238 

TOTAL a DATA POINTS - 2531 * DATA POINTS THIS UNLOADING - 147 

CORRELATION - .999677856393 

COD COMPLIANCE * 3.20855793284E-5 mm/N 

AREA - 9678.11102926 N»mm 
J - 198.876552682 KJ/m~2 

MAX LOAD LINE DEFLECTION - 1 . 26297067275 mm 
CHANGE IN COD - .0004572 mm 



UNLOADING * 8 

CRACK LENGTH - 15.0148445235 mm DELTA _a - .582163973634 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.0416778385 TO 14.9881516089 

TOTAL * DATA POINTS » 2761 * DATA POINTS THIS UNLOADING - 97 

CORRELATION - .999364204197 

COD COMPLIANCE - 3 .34673392487E-5 mm/N 

AREA * 10586.6765231 N«mm 
J = 215.448776655 KJ/m'2 

MAX lCAD LINE DEFLECTION - 1.36396114491 mm 
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CHANGE IN COD - .0006604 mm 



UNLOADING * 9 

CRACK LENGTH - 15.1216765395 mm DELTA _a - .689000984704 mm 

95% CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15 . 13586061 7 1 TO 15.1075318754 

TOTAL 0 DATA POINTS - 3042 * DATA POINTS THIS UNLOADING - 119 

CORRELATION - .99977636544 

COD COMPLIANCE - 3. 433430 78736E-5 mm/N 

AREA - 11709.4655627 N*mm 
J - 237.699077286 KJ/m*2 

MAX LOAD LINE DEFLECTION - 1.49020183752 mm 
CHANGE IN COD - .0005842 mm ' 



UNLOADING * 10 

CRACK LENGTH « 15.2624877891 mm DELTA _a > .829812234314 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.2760292626 TO 15.2489823604 

TOTAL 0 DATA POINTS - 3394 0 DATA POINTS THIS UNLOADING - 167 

CORRELATION - .999700915055 

COD COMPLIANCE - 3 .56785340653E-5 mm/N 

AREA - 12983.4083914 N-mm 
J * 262.394826944 KJ/m‘2 

MAX LOAD LINE DEFLECTION - 1.63760345064 mm 
CHANGE IN COD - .0007366 mm 



UNLOADING #11 

CRACK LENGTH * 15.3993548961 mm DELTA _a - .966679341294 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.4107739959 TO 15.3879615278 

TOTAL » DATA POINTS - 3867 * DATA POINTS THIS UNLOADING - 204 

CORRELATION - .999734089719 

COD COMPLIANCE - 3.70351235759E-5 mm/N 

AREA » 14557.9277896 N»mm 
J - 294.038000741 KJ/m‘2 

MAX LOAD LINE DEFLECTION - 1.82262379631 mm 
CHANGE IN COD - .0004318 mm 



UNLOADING * 12 

CRACK LENGTH - 15.669999544 mm DELTA _a - 1.2373239892 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.5847004127 TO 15.6553416287 

TOTAL 0 DATA POINTS - 4267 0 DATA POINTS THIS UNLOADING - 154 

CORRELATION - .999652892868 

COD COMPLIANCE - 4.00050024785E-5 mm/M 

AREA = 16330.4160222 N-mm 
J - 325.442585526 KJ/m'2 

MAX llIAD lINE DEFLECTION - 2.03825076401 nn 
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CHANGE IN COD - .0013462 mm 



UNLOADING » 13 

CRACK LENGTH - 15.8781418346 mm DELTA _a = 1.44546627973 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 15.8901071812 TO 15.8662051558 

TOTAL * DATA POINTS =4710 * DATA POINTS THIS UNLOADING = 208 

CORRELATION = .999677232181 

COD COMPLIANCE » 4.24377123766E-5 mm/N 

AREA = 17933.7183158 N-mm 
J » 355.360006494 KJ/m~2 

MAX LOAD LINE DEFLECTION = 2.24213294335 mm 
CHANGE IN COD » .000508 mm 



UNLOADING * 14 

CRACK LENGTH = 16.1357477216 mm DELTA _a = 1.7030721668 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.1463233402 TO 16.1251947229 

TOTAL * DATA POINTS = 5206 a DATA POINTS THIS UNLOADING = 239 

CORRELATION = .999696024845 

COD COMPLIANCE - 4 .57508725867E-5 mm/N 

AREA » 19646.0705407 N»mm 
J = 385.620902249 KJ/m‘2 

MAX LOAD LINE DEFLECTION = 2.47449867471 mm 
CHANGE IN COD = .0004826 mm 



UNLOADING a 15 

CRACK LENGTH = 16.4193583405 mm DELTA _a = 1.98668278568 mm 

957 CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.4364102869 TO 16.4023658462 

TOTAL a DATA POINTS » 5648 a DATA POINTS THIS UNLOADING = 194 

CORRELATION = .999324290353 

COD COMPLIANCE = 4.99261698493E-5 mm/N 

AREA » 21264.009056 N-mrn 
J = 412.083508069 KJ/m~2 

MAX LOAD LINE DEFLECTION = 2.70365472703 mm 
CHANGE IN COD » .0005334 mm 



UNLOADING a 16 

CRACK LENGTH = 16.9028626935 mm DELTA _a = 2.47018713865 mm 

957. CONFIDENCE INTERVAL ON CRACK LENGTH IS FROM 16.9171 1 T3874 TO 16.388656585 

TOTAL a DATA POINTS = 6102 a DATA POINTS THIS UNLOADING • 162 

CORRELATION • .999564028818 

COD COMPLIANCE = 5 . 8035G094828E-5 mm/N 

AREA - 22950.1834392 N*mm 
J = 428.162914109 Kj/m"2 

MAX LOAD LINE DEFECTION • 2.96634800605 mm 
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DATA FROM J vs DELTA.a PLOT 

SLOPE » 164.148013927 Y-INTERCEPT 
CORRELATION - .99667195294 FROM 9 



J 1 c 


a 


142.05774949 




TEST IS OVER! 




UNLOADING 


2 


DELTA a(mm> 


J(KJ/m 


1 




.016075 


14.8 


2 




.038476 


35.3 


3 




.054841 


63.3 


4 




.108820 


109.0 


5 




.308074 


169.9 


6 




• .382583 


185.3 


7 




.436845 


198:9 


8 




.532169 


215.4 


9 




.689001 


237.7 


10 




.829812 


262.4 


11 




.966679 


294.0 


12 




1 .237324 


325.4 


13 




1.445456 


355.4 


l-i 




1 .703072 


385.6 


15 




1.986683 


412.1 


16 




2.470187 


428.2 


DATA FROM 


J 


vs DELTA_a PLOT 




SLOPE -164, 


,148013927 Y- 


INTERCEPT 


CORRELATION 


- .99667195294 


FROM 9 



J 1 c - 142.05774949 



TEST IS OVER! 



- 123.933075018 
DATA POINTS. 



2 ) 



- 123.933075018 
DATA POINTS. 
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APPENDIX I 



J_DATA OUTPUT 

NHAT DATA FILE DO YOU WISH TO READ? 

J-INTEGRAL TEST OF HY80-5B BY: KURT 
DATE: 12-13-83 TIME: 10:35 

aO- 14.4326755548 WITH 95* CONDIFENCE INTERVAL 
FROM 14.4442813846 TO 14.42109582 

FROM COMPLIANCE = 2 . 84829954475E-5 

CORRELATION- .99984474937 FROM 176 DATA POINTS 

DO YOU WANT A HARD COPY OF THIS INFORMATION? 



FROM UNLOADING * 1 RECORDED 128 DATA POINTS 

COD COMPLIANCE - 2 .8755529381 7E-5 mm/N CORRELATION - .999529434569 
CRACK LENGTH - 14.4487507815 mm DELTA_A - .016075226716 mm 
95* CONFIDENCE INTERVAL FROM 14.4693514291 mm TO 14.4282322251 mm 
AREA - 711.365042717 N-mm J - 14.8049806501 KPa*m 
DATA POINT INDEXES FROM 504 TO 632 WERE USED. 



FROM UNLOADING * 2 RECORDED 132 DATA POINTS 

COD COMPLIANCE = 2.89307913363E-5 mm/N CORRELATION = .999467111958 
CRACK LENGTH - 14.4711513219 mm DELTA_A = .0384757670308 mm 
957. CONFIDENCE INTERVAL FROM 14.4927078135 mm TO 14.4496847244 mm 
AREA » 1693.69098454 N-mm J - 35.2970225265 KPa-m 
DATA POINT INDEXES FROM 815 TO 947 WERE USED. 



FROM UNLOADING * 3 RECORDED 157 DATA POINTS 

COD COMPLIANCE - 2.90942709989E-5 mm/N CORRELATION - .999017363069 
CRACK LENGTH = 14.4875164202 mm DELTA_A - .0548408654168 mm 
95* CONFIDENCE INTERVAL FROM 14.5143346851 mm TO 14.460837178 mm 
AREA - 3033.0003923 N-mm J - 63.28381022 KPa«m 
DATA POINT INDEXES FROM 1143 TO 1300 WERE USED. 



FROM UNLOADING * 4 RECORDED 109 DATA POINTS 

COD COMPLIANCE - 2.944429881 18E-5 mm/N CORRELATION - .999675611283 
CRACK LENGTH - 14.5414951598 mm DELTA_A - .108819604992 mm 
95* CONFIDENCE INTERVAL FROM 14.5601316376 mm TO 14.5229259715 mm 
AREA = 5208.32927275 N-mm J - 108.998439815 KPa-m 
DATA POINT INDEXES FROM 1561 TO 1670 WERE USED. 



FROM UNLOADING * 5 RECORDED 95 DATA POINTS 

COD COMPLIANCE - 3. 1 0248861 726E-5 mm/N CORRELATION = .999731 1 05207 
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CRACK LENGTH - 14.7407494376 mm DELTA_A - .308073882773 mm 
957. CONFIDENCE INTERVAL FROM 14.7587196422 mm TO 14.7228419911 mm 
AREA = 3216.7400787 N-mm J * 169.858689852 KPa-m 
DATA POINT INDEXES FROM 1965 TO 2060 WERE USED. 



FROM UNLOADING * 6 RECORDED 138 DATA POINTS 

COD COMPLIANCE = 3. 1 771 921 9461 E-5 mm/N CORRELATION = .998642018089 
CRACK LENGTH = 14.8152586049 mm DELTA_A * .382583050059 mm 
957. CONFIDENCE INTERVAL FROM 14.8481840789 mm TO 14.782543466 mm 
AREA = 9001.97373669 N-mm J = 135.342172083 KPa-m 
DATA POINT INDEXES FROM 2142 TO 2280 WERE USED. 



FROM UNLOADING * 7 RECORDED 147 DATA POINTS 

COD COMPLIANCE * 3 .20855798284E-5 mm/N CORRELATION = .999677856393 
CRACK LENGTH = 14.3695206824 mm DELTA_A > .436845127588 mm 
957 CONFIDENCE INTERVAL FROM 14.8849566722 mm TO 14.8541311238 mm 
AREA - 9678.11102926 N-mm J * 198.876552682 KPa-m 
DATA POINT INDEXES FROM 2384 TO 2531 WERE USED. 



FROM UNLOADING * 8 RECORDED 97 DATA POINTS 

COD COMPLIANCE - 3 . 34673392487E-5 mm/N CORRELATION » .999364204197 
CRACK LENGTH = 15.0148445285 mm DELTA.A - .582168973634 mm 
957 CONFIDENCE INTERVAL FROM 15.0416778385 mm TO 14.9881516089 mm 
AREA - 10586.6765231 N-mm J - 215.448776655 KPa-m 
DATA POINT INDEXES FROM 2664 TO 2761 WERE USED. 



FROM UNLOADING * 9 RECORDED 119 DATA POINTS 

COD COMPLIANCE » 3 .43343078736E-5 mm/N CORRELATION - .99977636544 
CRACK LENGTH » 15.1216765395 mm DELTA_A - .689000984704 mm 
957 CONFIDENCE INTERVAL FROM 15.1358606171 mm TO 15.1075318754 mm 
AREA - 11709.4655627 N-mm J - 237.699077286 KPa-m 
DATA POINT INDEXES FROM 2923 TO 3042 WERE USED. 



FROM UNLOADING * 10 RECORDED 167 DATA POINTS 

COD COMPLIANCE = 3 .56785340653E-5 mm/N CORRELATION * .999700915055 
CRACK LENGTH - 15.2624877891 mm DELTA_A = .829812234314 mm 
957 CONFIDENCE INTERVAL FROM 15.2760292626 mm TO 15.2489823604 ram 
AREA * 12983.4083914 N-mm J ■ 262.394826944 KPa-m 
DATA POINT INDEXES FROM 3227 TO 3394 WERE USED. 



FROM UNLOADING * 11 RECORDED 204 DATA POINTS 

COD COMPLIANCE * 3 . 70351235759E-5 mm/N CORRELATION = .999734089719 
CRACK LENGTH = 15.3993548961 mm DELTA_A - .966679341294 mm 
957 CONFIDENCE INTERVAL FROM 15.4107739959 mm TO 15.3879615278 mm 
AREA « 14557.9277896 N-mm J = 294.038000741 KPa-m 
DATA POINT INDEXES FROM 3663 TO 3867 WERE USED. 



FROM UNLOADING # 12 RECORDED 154 DATA POINTS 

COD COMPLIANCE - 4 . 00050024785E-5 mm/N CORRELATION = .999652892868 
CRACK LENGTH - 15.669999544 mm DELTA_A = 1.2373239892 mm 
957 CONFIDENCE INTERVAL FROM 15.6847004127 mm TO 15.6553416287 mm 
AREA - 16330.4160222 N-mm J = 325.442585526 KPa-m 
DATA POINT INDEXES FROM 4113 TO 4267 WERE USED. 
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FROM UNLOADING * !3 RECORDED 208 DATA POINTS 

COD COMPLIANCE - 4 .243771 2376SE-5 mm/N CORRELATION = .999677232181 
CRACK LENGTH = 15.8781418346 mm DELTA_A = 1.44546627973 mm 
957 CONFIDENCE INTERVAL FROM 15.8901071812 mm TO 15.8662051558 mm 
AREA - 17933.7183158 N*mm J - 355.360006494 KPa*m 
DATA POINT INDEXES FROM 4502 TO 4710 WERE USED. 



FROM UNLOADING * 14 RECORDED 239 DATA POINTS 

COD COMPLIANCE - 4 .57508725867E-5 mm/N CORRELATION - .999696024845 
CRACK LENGTH - 16.1357477216 mm DELTA A - 1.7030721668 mm 
957. CONFIDENCE INTERVAL FROM 16.1463233402 mm TO 16.1251947229 mm 
AREA - 19646.0705407 N*mm J - 385.620902249 KPa*m 
DATA POINT INDEXES FROM 4967 TO 5206 WERE USED. 



FROM UNLOADING * 15 RECORDED 194 DATA POINTS 

COD COMPLIANCE = 4 .99261698493E-5 mm/N CORRELATION - .999324290353 
CRACK LENGTH = 16.4193583405 mm DELTA_A - 1.98668278568 mm 
957 CONFIDENCE INTERVAL FROM 16.4364102869 mm TO 16.4023658462 mm 
AREA - 21264.009056 N»mm J - 412.083508069 KPa*m 
DATA POINT INDEXES FROM 5454 TO 5648 WERE USED. 



FROM UNLOADING * 16 RECORDED 162 DATA POINTS 

COD COMPLIANCE - 5 .80350094828E-5 mm/N CORRELATION « .999564028818 
CRACK LENGTH - 16.9028626935 mm DELTA_A = 2.47018713865 mm 
957 CONFIDENCE INTERVAL FROM 1 6.91 71 1 13874 mm TO 16.888656585 mm 
AREA - 22950.1834392 N«mm J = 428.162914109 KPa*m 
DATA POINT INDEXES FROM 5940 TO 6102 WERE USED. 
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DATA FROM aU rGR:HY80-5B 3Y: :<URT 
TESTED ON: 12-13-83 TIME: 10:35 



P r * 


LOAD < N ) 


DEFLECTION (mm) 


COD <mm) 


1 


2262.5 


. 079424 


.060576 


2 


2261.7 


.074706 


.060576 


3 


2265.2 


.076740 


.060637 


4 


2277.3 


.079003 


.061003 


5 


2350.8 


.093469 


.063539 


6 


2409.0 


.094973 


.065014 


7 


2442.3 


.095827 


.065395 


8 


2463.4 


.096491 


.066429 


9 


2478.9 


.096963 


.066954 


10 


2496.2 


.097378 


.067230 


1 1 


2505.2 


.097666 


.067627 


12 


2515.4 


.098165 


.067772 


13 


2526.1 


.093554 


.063207 


14 


2539.6 


.099114 


.063511 


15 


2554.0 


.099864 


.068969 


16 


2570.4 


. 100177 


.06941 1 


17 


2595.3 


.101411 


.070231 


18 


2624.1 


.102889 


.071 107 


19 


2653.4 


.103775 


.072045 


20 


2688.5 


. 104673 


.072903 


21 


2722.7 


.106054 


.073863 


22 


2746.9 


.106721 


.074519 


23 


2764.1 


.107052 


.075024 


24 


2779.2 


. 108077 


.075443 


25 


2799.9 


.108725 


.075956 


26 


2820.2 


.109596 


.076655 


27 


2837.3 


.110482 


.077142 


28 


2854.0 


.110667 


.077617 


29 


2873.4 


.112035 


.073313 


30 


2939.6 


. 114805 


.080284 


31 


2984.3 


.115696 


.081471 



l { l f 



6085 


5686.4 


2.954239 


2.460422 


6086 


5692.7 


2.954786 


2.460701 


6087 


5699.1 


2.955204 


2.461209 


6038 


5704.9 


2.955629 


2.461616 


6089 


5716.5 


2.956736 


2.462454 


6090 


5731 .8 


2.957790 


2.463294 


6091 


5747.5 


2.959093 


2.464359 


6092 


5761 .2 


2.959789 


2.465172 


6093 


5768.6 


2.960445 


2.465730 


6094 


5777.0 


2.961089 


2.466442 


6095 


5785.5 


2.961730 


2.467051 


6096 


5792.8 


2.962262 


2.467534 


6097 


5803.0 


2.963643 


2.463321 


6098 


5812.6 


2.964142 


2.469007 


6099 


5818.6 


2.964692 


2.469540 


6100 


5824.1 


2.965503 


2.470175 


6101 


5830.0 


2.965801 


2.470582 


6102 


5836.3 


2.966848 


2.471293 
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